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ABSTRACT 

This  thesis  explores  the  feasibility  of  making  wind  speed  profile  measurements  from  an 
oceanographic  buoy  using  a  Doppler  sodar. 

In  the  fall  of  2005,  we  deployed  a  Scintec  SFAS  sodar  on  an  ASIS  buoy.  Roughly  one 
week  of  buoy  motion  data  and  one  day  of  sodar  observations  were  collected.  Data  from 
both  this  deployment,  and  the  Martha’s  Vineyard  Coastal  Observatory,  were  used  in 
conjunction  with  models  to  predict  sodar  performance.  Results  are  compared  for  an  ASIS 
and  a  3-meter  discus  buoy.  We  also  predict  the  yearly  average  probability  of  sodar  data 
availability  in  the  presence  of  buoy  motion.  We  show  that  buoy  tilting  in  response  to 
wave  forcing  is  the  main  factor  affecting  sodar  performance.  Our  results  strongly  suggest 
that  ASIS  is  a  suitable  platform  for  sodar  measurements  at  sea. 


Thesis  Supervisor:  Dr.  Eugene  A.  Terray 
Title:  Research  Specialist 
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Chapter  1 
Problem  Overview 
1.1  Motivations 

In  the  present  day,  fixed  meteorological  towers  are  the  most  commonly  used  platfonn 
from  which  to  make  wind  measurements  in  the  marine  atmospheric  boundary  layer 
(MABL).  The  determining  factors  in  the  cost  of  a  meteorological  tower  are  the  depth  of 
the  water,  the  air-side  height  of  the  tower  and  the  typical  wave  climate  of  the  chosen 
location.  As  a  result,  it  is  not  economical  to  build  a  tower  that  extends  far  above  the  mean 
water  level  or  one  that  is  situated  in  deep  water.  For  example,  it  can  cost  between  $1.5  - 
2  M  to  construct  a  60  m  tower  in  10  m  of  water.  From  the  point  of  view  of  site 
characterization  for  wind  power  generation,  fixed  meteorological  towers  clearly  are 
limited  to  measuring  winds  at  a  fixed  location,  so  there  is  no  possibility  of  exploring 
different  or  better  sites.  Also,  towers  generally  are  not  high  enough  to  measure  winds 
over  the  span  of  large  wind  turbine  rotors.  For  example,  the  General  Electric  3.6sl 
offshore  wind  turbine  has  a  rotor  diameter  of  1 1 1  m.  Fixed  meteorological  towers  also 
are  inadequate  when  it  comes  to  capturing  the  full  scope  of  MABL  physics.  In  the 
summer  months  along  the  Atlantic  continental  shelf,  stable  stratification  in  the 
atmosphere  can  lead  to  the  formation  of  low-level  jets  (Vickers  and  Mahrt,  2004).  An 
example  is  shown  in  Figure  1.1,  which  is  taken  from  the  CBLAST  -  Low  experiment 
conducted  off  the  south  coast  of  Martha’s  Vineyard,  Massachusetts  in  August  2001.  The 
wind  speed  profile  shows  a  strong  shear  over  the  lower  100  m,  reaching  a  maximum  of 
about  14  m/s  at  1 10  m,  where  it  reverses  direction.  Because  such  jets  can  produce  loads 
that  will  shorten  the  life  of  wind  turbine  blades,  it  is  important  to  be  able  to  characterize 
them  and  determine  their  frequency  of  occurrence,  and  hence  wind  profiles  over  heights 
of  order  150  m  are  necessary. 

In  this  thesis  we  will  study  the  feasibility  of  using  Doppler  sodar,  deployed  on  an 
oceanographic  buoy,  to  profile  winds  over  several  hundred  meters  with  a  range  resolution 
of  several  meters. 
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Figure  1.1  Data  taken  from  the  CBLAST-  Low  experiment  conducted  off  the  south 
coast  of  Martha’s  Vineyard  in  August  2001  (Vickers  and  Mahrt,  2004). 

1.2  Doppler  profiling  instruments 

In  this  section  we  survey  approaches  to  measuring  wind  profiles  over  several  hundred 
meters.  Three  different  kinds  of  instruments  are  capable  in  principle  of  sensing  wind 
profiles:  Doppler  lidars,  radars,  and  sodars. 

Doppler  lidar  (laser  radar)  measures  the  Doppler  shift  of  radiation  that  is  scattered  by 
natural  particles  such  as  dust  and  water  droplets. 

Two  primary  varieties  of  lidars  exist:  pulsed  and  continuous  wave  (CW).  Pulsed  lidars 
have  an  along  -  beam  range  resolution  that  is  constant  with  respect  to  range,  but  typically 
cost  $1M  or  more  because  the  pulsed  lasers  and  digitizer  are  expensive.  In  contrast,  CW 
lidars  tend  to  be  much  less  expensive.  For  example,  QinetiQ  makes  a  CW  infrared  lidar 
called  the  ZephIR  (see  Figure  1.2),  that  is  based  on  readily  available  parts  from  the 
telecom  industry,  and  sells  for  ~$180K.  It  uses  a  conically  scanned  laser  at  a  wavelength 
of  1.55  pin.  The  ZephIR  is  a  compact  package  that  weighs  40  kg  and  consumes  about 
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100  W.  The  system  can  make  wind  speed  measurements  up  to  200  m,  and  data  from  it 
has  compared  very  well  with  cup  anemometer  measurements  (Jorgenson  et.  al.,  2004). 
The  main  difficulty  in  the  CW  approach  is  how  to  determine  range.  The  ZephIR  uses  a 
mechanically  scanned  focal  point.  However,  in  this  case  the  range  cell  increases  as  the 
square  of  the  range.  So,  for  example,  whereas  at  50  m  the  range  resolution  is  4  m,  at 
200  m  the  range  resolution  has  decreased  to  60  m  (Jorgenson  et.  al.,  2004). 


Figure  1.2  ZephIR  lidar  and  its  principle  of  operation  (figure  courtesy  of  D.A. 

Smith,  QinetiQ  Ltd.) 


A  second  way  to  obtain  wind  profiles  in  the  MABL  is  via  Doppler  radar.  These  devices 
emit  electromagnetic  energy  into  the  atmosphere  and  detect  the  Doppler  shift  of  the 
returned  signal.  NOAA  currently  operates  a  network  of  Doppler  radars  in  the  United 
States  that  profile  winds  up  to  several  kilometers.  Although  these  instruments  are 
commercially  available,  they  have  a  range  resolution  of  60-100  m,  making  them 
unsuitable  for  our  application.  At  the  present  time,  there  are  no  commercially-available 
radars  having  a  range  resolution  of  order  a  few  meters  over  the  lowest  200  m. 

ProSensing  has  developed  a  prototype  radar  called  the  MBL-profiler  (see  Figure  1.3)  that 
can  be  used  from  a  small  boat  or  buoy.  The  MBL-profiler  uses  a  four-beam  microstrip 
antenna  and  can  probably  profile  over  several  hundred  meters.  Although  not 
commercially  available,  when  in  regular  production,  such  a  radar  would  likely  cost 
O($200K). 
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Figure  1.3  ProSensing's  MBL-profiler  rigged  aboard  a  small  boat 

Doppler  sodars  work  by  scattering  an  acoustic  pulse  from  sound  speed  fluctuations  in  the 
atmosphere  and  estimating  wind  speed  at  each  range.  Wind  profiles  are  computed  from 
the  Doppler  shift  which  is  computed  from  the  frequency  of  the  returned  signal  (more 
specifics  are  discussed  later).  Several  compact,  low-power,  sodars  are  available 
commercially.  For  example,  Scintec  and  Remtech  both  build  flat  phased  array  sodars  that 
require  — 100  W  (transmit  power).  These  instruments  can  profile  over  ~  150  m  with  5  m 
range  resolution  and  are  comparatively  priced  at  roughly  $  60  K.  Since  we  have  access 
to  the  Scintec  SFAS  model  we  will  focus  on  that  sodar  here. 
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Chapter  2 
Sodars 

2.1  Sodar  research  background 

Sodars  have  a  long  history  of  use  on  land  and  have  been  successfully  used  in  many 
pollution  and  boundary  layer  studies.  Numerous  comparison  experiments  have  been 
performed  to  validate  the  accuracy  of  Doppler  sodar  use  over  land  (Crescenti,  1997). 
However,  few  marine-based  sodar  studies  have  been  conducted  and  the  results  have  been 
mixed.  It  has  been  suggested  that  one  reason  for  the  poor  results  reported  is  that  the 
marine  environment  may  not  provide  enough  scattering  strength  (Fairall  et.al.,  1997). 
However,  in  a  1  day  study,  Barthehnie  et.  al.  (2003)  mounted  a  sodar  on  a  moored  ship 
near  Vindeby,  Denmark  and  got  good  results.  Their  Aerovironment  4000  mini-sodar 
obtained  results  that  were  not  significantly  different  than  those  collected  during  land- 
based  experiments  with  the  same  instrument  (Coelingh  et.  al.  2000). 

A.  L.  Rogers  (2006)  has  also  obtained  good  quality  results  using  a  sodar  in  the  marine 
environment.  In  the  fall  of  2005,  his  group  operated  their  ART-VT1  sodar  from  the 
WHOI  Air  Sea  Interaction  Tower  (ASIT)  south  of  Martha’s  Vineyard.  The  sodar  was 
deployed  from  22  September  2005  to  18  November  2005  and  got  good  data  return. 

Figure  2.1  shows  a  plot  of  percentage  of  acceptable  10-minute  averages  as  a  function  of 
range.  For  heights  up  to  80  m  above  the  sodar,  at  least  50%  of  the  10  minute  averaged 
spectra  yielded  acceptable  Doppler  estimates. 

The  experiences  of  Barthehnie  et.  al.  (2003)  and  A.  L.  Rogers  (2006)  suggest  that  there  is 
potential  for  obtaining  quality  data  from  a  sodar  in  the  marine  enviromnent.  However,  the 
sodar  studies  performed  by  the  previously  mentioned  groups  were  both  conducted  from 
platforms  where  motion  was  negligible.  Rogers  (2006)  recorded  sodar  data  from  a  tower, 
whereas  Barthehnie  et.  al.  (2003)  made  measurements  from  a  ship  anchored  with  a  three- 
point  mooring  in  low  sea  states.  No  one  has  attempted  to  collect  sodar  data  from  a 
moving  platfonn  in  the  marine  environment.  This  thesis  explores  the  feasibility  of 
obtaining  quality  wind  profile  data  from  an  oceanographic  buoy. 
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Figure  2.1  Percentage  of  acceptable  10  minute  averages  (A.  L  Rogers,  2006). 


2.2  Sodar  theory  of  operation 

A  sodar  works  by  emitting  an  acoustic  pulse  that  is  backscattered  from  sound  speed 
inhomogeneities  in  the  atmosphere.  Monostatic  sodars  have  their  transmitter  and  receiver 
in  the  same  location  so  that  only  direct  backscatter  can  be  measured.  Backscatter  is 
detected  if  turbulent  eddies  fulfill  the  Bragg  condition  ( i.e .  the  spatial  wavenumber  of  the 
turbulent  temperature  fluctuations  twice  the  acoustic  wavenumber).  Sound  can  also  be 
backscattered  from  density  discontinuities.  The  Doppler  shift  is  computed  from  the  return 
frequency  of  the  backscattered  signal,  The  magnitude  of  the  frequency  shift  is 


=  -2— /e  (2.1) 
c 


where  fe  is  the  frequency  of  the  emitted  pulse,  ur  is  the  component  of  the  wind  velocity 
along  the  beam,  and  c  is  the  speed  of  sound  (Pierce,  1981).  For  a  beam  tilted  at  an  angle 
0  with  respect  to  the  vertical,  the  Doppler  shift  is 


Acceptable  10-minute  Averages 
September  22  -  November  18,  2005 
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(2.2) 


Of  = -2  — fe  sin(0)  -2  — fe  cos(0) 
c  c 


where  w  and  w  are  the  horizontal  and  vertical  components  of  the  wind  in  the  plane 
defined  by  the  vertical  and  the  beam.  Using  (2.2),  easterly  (u),  northerly  (v),  and  vertical 
wind  components  (w)  can  be  calculated  from  frequency  shifts  measured  by  north  -  (N), 
east  -  (E)  and  vertical-pointing  beams  (V)  as 


u  (>f  \:  c  ,  ¥y  C 

2 fE  sm(6E)  2 fv  tan(0£) 


(2.3) 


C  t  Wy  c 

2 fN  sin(0w)  2 fv  tan(0,v) 


(2.4) 


w  =  - 


2  fv 


(2.5) 


2.3  Spectral  first  moment 

This  section  describes  how  to  estimate  Doppler  shift  numerically  via  a  spectral  first 
moment  calculation.  With  an  estimate  of  the  spectral  mean  frequency  (/),  we  compute 
Doppler  shift  using  the  following 


$ f-f'-f- 


(2.6) 


The  first  moment  of  the  measured  spectrum,  SM  ,  yields  the  mean  frequency,  j) 


Ib  = 


f  fSM(f)df  f(f  2  fN^SMdf  l 


/  SM(f)df  f  SMdf 


+  2  f" 


(2.7) 
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where  SM  is  the  measured  spectrum  of  the  returned  pulse  and  fN  is  the  Nyquist 
frequency.  This  estimate  is  biased  because  the  spectrum  of  the  noise  is  included  in  the 
calculation. 


Equation  (2.7)  can  be  written  as 


fa  =«/  +  ^(1-«)/jv 

where 

a.l—L.  1  JjL,  1  -L _ SNR» -1 

fsMdf  JV„  SNR,  SNR„ 
From  equation  (2.8),  we  estimate  the  mean  frequency  as 

f  =  UB  -^(l-a)/jv]/a. 


(2.8) 


(2.9) 


(2.10) 


2.4  Variance  in  the  Doppler  frequency  estimate 

In  this  section,  we  use  the  ratio  of  signal  over  noise  to  assess  the  quality  of  our  sodar  data. 
For  a  given  SNR  we  can  calculate  how  much  variance  we  have  in  the  Doppler  frequency 
estimate.  The  Cramer-Rao  bound  for  an  unbiased  estimator  is: 

of  =(2jtT)-\l+ 36/ SNR +  30 /SNR2)112  (2.11) 

(Theriault,  1986)  where  of  is  Doppler  frequency  variance,  T  is  the  pulse  length  and  SNR 

is  the  signal-to-noise  ratio.  The  variance  of  Doppler  frequency,  non-dimensionalized  by 
the  pulse  length,  T,  and  multiplied  by  2n,  is  shown  in  Figure  2.2.  As  SNR  increases, 
OjlnT  decreases  rapidly  and  asymptotes  to  unity. 
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Non-dimensionalized  Doppler  frequency  variance 


Figure  2.2  Non-dimensionalized  Doppler  frequency  variance  versus  SNR.  The 

range  resolution  is  A R  =  ^  CT . 

ForC  =  330m/s,  a30ms  pulse  has  a  range  resolution  of  5  m.  Hence  for  SNR  between 
one  and  infinity,  of  ~  5  ~  40 Hz .  The  corresponding  uncertainty  for  fe~  330  Hz  in  the 

radial  velocity  is  0.25  -  2  m/s. 

2.5  Scintec  SFAS 

The  Scintec  SFAS  is  a  compact  sodar  capable  of  profiling  winds  up  to  200  m  with  0(5m) 
range  resolution.  The  sodar,  shown  in  Figure  2.3,  is  a  64  element  array  that  can 
synthesize  up  to  nine  shaded  beams  in  transmit  and  receive  mode  at  angles  of  0°,  ±19° 
and  ±24° .  The  SFAS  can  emit  up  to  10  tones  (acoustic  pulses)  in  the  frequency  range  of 
2525  to  4850  Hz.  The  instrument  is  relatively  portable  in  that  it  weighs  only  1 1.5  kg,  and 
measures  44  x  42  x  16  cm.  It  sits  in  a  baffle  (Figure  B.2)  that  opens  at  an  angle  of  25.9°. 
The  baffle  serves  to  mitigate  the  effect  of  high  angle  noise  on  the  received  signal. 
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Figure  2.3  Scintec  SFAS  sodar. 


2.6  Scintec  beampattern  calculation 

In  this  section  we  compute  the  Scintec  SFAS  beampattern.  A  description  of  the  relevant 
theory  is  given  in  Appendix  A.  The  directivity  of  the  vertical  (0°)  and  steered  beams  (± 
19°,  ±  24°)  is  slightly  different.  We  compare  a  1 -dimensional  slice  through  the 
beampatterns  for  the  vertical  and  steered  19°  beam  case  at  the  instrument’s  middle  tone, 
3568  Hz  (see  Figure  2.4).  In  the  figure,  the  steered  beampattern  for  the  19°case  is 
centered  at  0°  so  that  it  can  easily  be  compared  to  the  vertical  beam.  Only  positive  angles 
of  the  beampattern  are  shown  as  the  beampattern  is  a  mirror  image  at  the  corresponding 
negative  angles.  Since  the  SFAS  beamforms  both  the  transmit  and  receive  beams,  we 
show  the  product  of  the  transmit  and  receive  beam  intensity  patterns  (i.e.  the  2-way 
intensity)  normalized  to  have  a  maximum  value  of  unity  along  the  central  axis. 
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Steered  and  Vertical  Unshaded  Beampatterns 


Figure  2.4  A  comparison  of  steered  and  vertical  beampatterns  for  the  unshaded 
case.  Beampatterns  are  plotted  as  normalized  two-way  intensity. 


The  vertical  and  steered  beampatterns  have  very  similar  main  lobes.  Table  2. 1 
summarizes  the  characteristics  of  the  unshaded  steered  and  vertical  beams  (the  19°  beam 
characteristics  are  given  relative  to  the  mainlobe). 


Beam 

19° 

0° 

Location  of  Sidelobes 

-71.4°, -19.6°, 19.6°, 71.4° 

-46°, -17. 2°, 17. 2°,  46° 

Location  of  Minima 

-39.3°, -1 1.9°,  11.9°,  39.3° 

-28°, -10.8°, 10.8°, 28° 

Relative  intensity  at  90° 

1.315  xl0“2 

5.373  xl0“3 

3dB  Width 

10° 

9.6° 

6dB  Width 

14° 

12.8° 

Table  2.1  Characteristics  of  unshaded  19°  and  0°  beampatterns. 


To  reduce  the  sidelobe  amplitude,  the  SFAS  can  also  beamform  using  shading.  A 
comparison  of  the  normalized  2-way  intensity  for  the  1 -dimensional  shaded  and 
unshaded  vertical  beam  is  shown  in  Figure  2.5.  Again,  only  positive  angles  are  shown. 
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Note  that  shading  decreases  the  intensity  of  the  sidelobes,  while  widening  the  main  lobe. 
It  is  also  noteworthy  that  the  shaded  beampattern’s  intensity  is  about  a  factor  of  5  higher 
than  the  unshaded  beampattern’s  intensity  at  90  degrees. 

1 D  Shaded  &  Unshaded  Intensities  for  f  =  3568  Hz 


Figure  2.5  Scintec  1-dimensional  beampattern  for  the  shaded  and  unshaded  vertical 
beampattern.  Intensity  is  two-way  and  normalized. 

The  important  characteristics  of  the  1 -dimensional  shaded  and  unshaded  vertical 
beampatterns  are  summarized  in  Table  2.2. 


Shaded  Beampattern 

Unshaded  Beampattern 

Location  of  Sidelobes 

-48°, -20°, 20°, 48° 

-46°, -17.2°, 17.2°, 46° 

Location  of  Minima 

-34°, -12.8°, 12.8°, 34° 

-28°, -10.8°, 10.8°, 28° 

Relative  intensity  at  90° 

0.0290 

5.373  xl0“3 

3dB  Width 

11.2° 

9.6° 

6dB  Width 

15.2° 

12.8° 

Table  2.2  Summary  of  Scintec  vertical  beampattern  characteristics. 
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The  normalized  two-way  intensity  as  a  function  of  ( x,y )  =  (kA /jt) si ndtcos cosine/))  for 
the  2-dimensional  shaded  SFAS  vertical  beampattern  is  shown  in  Figure  2.6.  Here  A  = 
0.05  m  is  the  transducer  spacing  in  both  x  and  y.  Note  that  the  central  lobe  is 
approximately  azimuthally  symmetric. 

2-d  Normalized  Intensities  at  f  =  3568  Hz 


-1  -0.5  0  0.5  1 

(  k  A  /  n  )  sin0  cosq 


Figure  2.6  2-dimensional  normalized  intensities  for  the  shaded  SFAS  vertical 
beampattern  at  3568  Hz  (k  ~  0.1m). 
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Chapter  3 
Sodar  Platforms 
3.1  Overview 

Since  buoy  motion  decreases  SNR  by  reducing  the  overlap  between  the  transmit  (T)  and 
receive  (R)  beams,  it  is  important  to  understand  how  the  shape  of  a  buoy  affects  its 
response  to  wave  forcing. 

We  want  to  choose  a  buoy  that  does  not  experience  a  significant  amount  of  rotational 
motion  (i.e.  pitch,  roll,  and  yaw).  Tilts  in  pitch  and  roll  will  likely  have  a  negative  effect 
on  sodar  performance  since  the  main  lobe  of  the  sodar’ s  beampattern  has  a  6-dB  width  of 
only  10°.  Hence  buoy  tilts  in  pitch  and  roll  need  to  be  small  compared  to  the  width  of  the 
main  lobe  in  order  for  the  sodar  to  get  a  strong  return.  This  is  a  stiff  requirement  since 
many  oceanographic  buoys  often  tilt  in  excess  of  10°  due  to  wave  forcing.  In  contrast, 
their  heading  changes  typically  are  slow,  so  that  rotational  motion  in  yaw  is  expected  to 
have  a  smaller  impact  on  sodar  performance  than  pitch  and  roll. 

Translational  motion  (heave,  surge,  and  sway)  will  also  have  a  negative  effect  on  sodar 
performance,  but  to  a  lesser  degree  than  rotational  motions.  Sodar  range  bins  are 
typically  5  -  10  m,  so  buoy  motions  in  heave,  for  example,  are  small  compared  to  sodar 
range  resolution.  Therefore,  there  is  not  as  strict  a  requirement  on  buoy  heave. 

At  one  end  of  the  spectrum  of  oceanographic  buoys  there  are  large  spar  buoys,  such  as 
the  110  m  long  “Floating  Instrument  Platform”  (FLIP),  operated  by  the  Scripps 
Institution  of  Oceanography,  which  does  not  respond  to  wave  forcing 
(http://sio.ucsd.edu/voyager/flip/index.html).  FLIP  is  stable  in  waves  because  its  resonant 
frequency  is  below  the  wave  band.  The  following  example  explains  why  spar  buoys  are 
stable.  A  cylinder  of  length,  L,  floating  upright  has  a  resonant  period  of  2jt^L/ g  . 

Hence,  a  simple  spar  greater  than  25  m  in  length  will  be  insensitive  to  10  s  waves,  and  a 
spar  buoy,  such  as  FLIP,  that  is  greater  than  100  m  long  will  show  no  response  to  20  s 
waves.  Spar  buoys  also  need  freeboard  so  that  they  are  not  overtopped  by  passing  waves. 
A  large  spar  would  be  a  nice  platform  to  host  a  sodar  if  one  is  looking  to  make 
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measurements  far  offshore.  From  the  point  of  view  of  wind  turbine  site  characterization 
and  coastal  atmospheric  physics,  however,  a  large  spar  buoy  will  not  fit  the  bill.  A  buoy 
that  is  not  portable  and  clearly  cannot  be  deployed  in  coastal  waters,  will  not  meet  our 
needs  since  we  require  a  platform  that  can  easily  be  deployed  in  a  variety  of  locations, 
including  littoral  areas. 

In  contrast  to  spar  buoys  is  the  discus  buoy,  which  follows  the  ocean  surface  in  pitch, 
roll,  and  heave.  These  buoys  also  tend  to  follow  the  surface  in  surge  and  sway, 
depending  on  the  mooring  arrangement.  An  example  of  a  discus  buoy  is  the  three-meter 
discus  buoy  used  by  the  National  Data  Buoy  Center  (NDBC)  to  measure  waves.  Since 
the  3 -meter  discus  buoy  is  so  widely  used  in  Oceanographic  research,  it  is  of  interest  to 
evaluate  its  suitability  as  a  possible  platform  in  comparison  to  other  platforms.  We  do 
this  in  Section  3.2,  below. 

A  buoy  design  that  shares  the  characteristics  of  a  stable  spar  buoy  as  well  as  a  discus 
buoy  is  a  hybrid  spar  buoy  called  the  “Air  Sea  Interaction  Spar”  (ASIS),  shown 
schematically  in  Figure  3.4.  ASIS  is  a  short  spar  with  a  heavily  damped  resonance.  As  a 
result,  its  response  to  waves  resembles  a  low-pass  filter  with  a  cutoff  frequency  ranging 
from  6  to  8  seconds,  depending  on  its  particular  configuration.  The  ASIS  buoy  is  part  of 
a  system  that  includes  a  50  m  floating  tether  that  connects  it  to  a  surface  mooring.  This 
allows  ASIS  to  follow  wave  motion  relatively  freely  in  heave,  surge  and  sway.  However, 
ASIS  tilts  in  pitch  and  roll  generally  are  small.  ASIS  is  potentially  a  good  platform  to 
host  a  sodar,  so  we  will  discuss  its  design  in  more  detail  in  Section  3.3. 

3.2  Three  meter  discus  buoy  characteristics 

Since  sodar  works  by  measuring  the  Doppler  shift,  the  mis  error  in  the  wind  velocity 
measurement  is  affected  by  the  motion  of  the  instrument.  In  this  section,  we  investigate 
the  rms  tilts  of  a  3-meter  discus  buoy.  A  discus  buoy  3  meters  long  follows  waves 
roughly  up  to  twice  its  diameter  (6  meters).  This  corresponds  to  a  wavenumber, 
k  =  1  m~l .  With  k  we  can  solve  the  dispersion  relation  for  deep  water,  a>  =  -y[gk ,  to  get 
the  maximum  resolvable  frequency  seen  by  the  buoy.  In  this  case,  a>  =  2 jtf=  3.3  rad/s,  so 
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a  3-m  discus  buoy  ideally  will  follow  the  surface  slope  of  waves  up  to  frequencies  (/)  of 
0.5  Hz. 

To  get  an  understanding  of  how  the  buoy  behaves  in  a  variety  of  different  conditions,  we 
look  at  probability  distributions  of  buoy  rms  tilts  for  2005.  To  predict  the  rms  tilts  of  the 
discus  buoy,  we  apply  a  transfer  function  to  a  year’s  worth  (2005)  of  Martha’s  Vineyard 
Coastal  Observatory  (MVCO)  wave  slope  spectra.  Our  transfer  function  has  an  ideal 
low-pass  filter  form  that  is  unity  from  0  to  0.5  Hz,  and  vanishes  at  higher  frequencies.  A 
plot  of  the  probability  distribution  of  rms  tilts  computed  over  20  minute  bursts  for  the 
year  2005  is  given  in  Figure  3. 1 


rms  tilts  computed  over  20  min  bursts  2005 


Figure  3.1  Distribution  of  3-meter  discus  buoy  rms  tilts  for  2005. 

Reductions  in  SNR  due  to  three  meter  discus  buoy  tilts  can  be  considerable,  so  we  want 
to  get  an  idea  of  how  other  buoys  respond  under  the  same  conditions.  To  do  this  we 
compare  the  performance  of  the  three  meter  discus  buoy  to  the  performance  of  an  Air-Sea 
Interaction  (ASIS)  Spar  for  the  same  days  of  data.  We  have  six  days  of  ASIS  motion 
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data  from  the  ASIS  that  we  deployed  south  of  the  MVCO  tower  in  the  fall  of  2005  (the 
details  of  this  experiment  are  described  in  Appendix  B,  while  the  calibrations  of  the 
sensors  used  in  the  experiment  can  be  found  in  Appendix  D).  In  figure  3.2,  we  plot  a 
comparison  of  measured  ASIS  rms  tilts  to  3 -meter  discus  buoy  tilts  (for  the  same  days 
that  we  have  ASIS  data),  predicted  from  MVCO  measurements.  On  average,  3-meter 
discus  buoy  rms  tilts  are  roughly  4  times  greater  than  ASIS  tilt.  Since  ASIS  buoy  tilts  are 
much  smaller,  we  continue  an  investigation  into  the  characteristics  of  the  ASIS  platform. 

Comparison  of  Buoy  rms  Tilts 


Figure  3.2  A  comparison  of  3-meter  discus  buoy  predicted  rms  tilts  to  ASIS 
measured  rms  tilts.  The  solid  line  shows  a  linear  least  squares  fit  to  the  data  and  has 

a  slope  of  3.8. 


3.3  ASIS  design  characteristics 

In  this  section,  we  discuss  the  general  design  principles  underlying  ASIS,  and  sketch  a 
typical  buoy  design  and  instrumentation  payload  to  support  sodar  measurements.  ASIS 
(Graber  et.al.,  2000)  is  a  short,  highly  damped  spar  buoy  (Figure  3.3).  Unlike  a  true  spar 
buoy,  ASIS  moves  with  long  waves  (z'.e.  periods  >  7  -  8  s)  but  is  relatively  insensitive  to 
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shorter  wavelengths.  Because  it  has  a  small  response  to  pitch  and  roll,  it  likely  is  a  good 
platform  for  air-sea  measurements.  ASIS  has  the  added  capability  of  being  able  to  be 
disassembled  and  packed  in  a  20-foot  shipping  container.  It  is  also  advantageous  in  that 
it  can  be  deployed  from  small  research  vessels. 


Data  Acquisition  System 
Temperature/humidity 
^Flotation 


Tether 


Radar  reflector 

Marine  light  |  SFAS  baffle  Radiometer 

-Scintec  SFAS  SODAR 
-ARGOS 


2.83m 


SFAS  processing  unit 


1.72m 


,24m 


Doppler  sonar 


4,34m 


Battery  cans  (5) 


Figure  3.3  Drawing  of  a  typical  ASIS  buoy. 

ASIS  is  designed  to  follow  long  waves  and  be  stable  in  short  waves.  Hence,  the  dynamic 
range  of  perceived  wave  motion  by  the  buoy  is  reduced,  thus  minimizing  the  risk  of 
overtopping  in  the  presence  of  long  waves  with  larger  amplitudes.  Because  the  buoy  has 
considerable  draft,  it  can  be  ballasted  to  be  stable  in  pitch  and  roll  in  typical  New 
England  sea  states.  The  buoy’s  buoyancy  is  primarily  provided  by  the  five  upright 
columns  and  the  five  floats  (horizontal  cylinders)  that  are  situated  at  the  base  of  the 
columns.  The  main  tube  is  free-flooding  and  houses  adjustable  ballast  for  trim.  ASIS  is 
typically  ballasted  so  that  the  mean  water  level  is  in  the  middle  of  the  five  uprights.  The 
foot  of  the  buoy  is  made  up  of  expanded  aluminum  grating  to  increase  drag  in  heave. 
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Battery  cans  are  located  at  the  foot  and  add  a  considerable  amount  of  ballast,  which  helps 
to  keep  the  buoy  stable. 


Although  ASIS  can  be  used  as  a  freely  drifting  surface  platform,  it  is  more  useful  in  our 
application  as  a  moored  buoy.  Since  the  five  uprights  do  not  provide  much  reserve 
buoyancy  (they  total  124  lbs/ft  of  immersion),  ASIS  must  be  tethered  using  a  floating  line 
to  a  nearby  surface  mooring  in  order  to  prevent  large  vertical  forces  on  the  buoy.  The 
ASIS  buoy,  tether,  and  surface  mooring  make  up  the  ASIS  “system”. 


The  dimensions  of  a  typical  ASIS  buoy  are  shown  in  Table  3.1.  The  corresponding  air 
weights  and  buoyancies  are  shown  in  Table  3.2.  The  draft  of  ASIS  can  adjusted  by 
adding  or  removing  weights  from  the  main  tube. 


Item 

Quantity 

Shape 

Dimensions  (LxD) 

(inches) 

Foot 

1 

Pentagon 

4.2x39 

Batteries 

5 

Cylinder 

38x12.7 

Main  Tube 

1 

Cylinder 

170x12 

Candelabra  tubes 

5 

Cylinder 

53x4 

Floats 

5 

Cylinder 

44x12.7 

Uprights 

5 

Cylinder 

137x8.4 

Top  Plate 

1 

Pentagon 

4.2x49 

Table  3.1  Dimensions  of  ASIS  principal  structural  components. 


The  dimensions  of  ASIS  can  be  varied  somewhat  in  order  to  fine-tune  its  response  to 
different  sea  states.  The  response  functions  of  ASIS  generally  have  the  shape  of  a  low 
pass  filter.  The  drag  due  to  the  expanded  metal  in  the  foot  damps  the  resonant  heave 
response  of  the  buoy  (Graber  et.  al.,  2000). 
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Item 

Details 

Quantity 

Weight  in  Air 

Buoyancy 

(pounds) 

(pounds) 

ASIS  Buoy 

Foot 

1 

118.9 

44.8 

Batteries 

5 

1115 

887.5 

Main  Tube 

1 

241.5 

91 

Candelabra  tubes 

5 

69 

123 

Floats 

5 

136 

1027 

Uprights 

5 

634.5 

700.5 

Top  Plate 

1 

266.1 

Ballast 

114.8 

Instrumentation 

SFAS  SODAR 

1 

46 

SFAS  table 

1 

36 

SFAS  processing  unit 

1 

38 

SFAS  baffle 

1 

20 

Data  Acquisition  Box 

1 

21 

ARGOS 

1 

5 

Freewave  modem 

1 

2 

Marine  Light 

1 

5 

Radar  Reflector 

1 

5 

Total  Weight 

2874 

2874 

Table  3.2  Weights  and  buoyancy  of  ASIS  structural  components  and  sensors. 


3.4  ASIS  instrumentation  and  data  acquisition 

In  addition  to  wind  profiling,  ASIS  is  capable  of  making  other  measurements.  These 
measurements  can  be  divided  into  three  categories:  meteorological,  buoy  motion,  and 
waves. 


Meteorological  measurements  of  interest  are  temperature,  relative  humidity,  solar 
radiation,  barometric  pressure  and  rainfall.  Eppley  Labs  and  RM  Young  manufacture 
many  of  the  instruments  capable  of  these  measurements.  A  typical  list  of  instrumentation 
is  given  in  Table  3.3.  Data  can  be  digitized  and  recorded  with  the  buoy’s  main  data 
acquisition  system  or  a  separate  data  logger  such  as  a  Campbell  CR-10X  logger. 

The  buoy’s  motion  is  measured  with  an  inertial  motion-sensing  unit  (IMU).  The  IMU 
consists  of  three  orthogonal  accelerometers  and  three  orthogonal  rate  gyros.  A  compass 
must  also  be  included  to  measure  buoy  heading.  BEI  Systron  Donner  makes  a  commonly 
used  IMU  (see  Table  3.3).  Since  accelerometers  cannot  distinguish  between  tilt  and 
lateral  acceleration,  rate  gyros  must  be  used  to  determine  tilt  in  the  band  of  wave 
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frequencies.  Rate  gyros  measure  angular  velocity  about  an  axis  and  must  be  integrated  to 
obtain  tilt.  At  frequencies  well  below  the  wave  band,  lateral  accelerations  are  negligible, 
so  accelerometers  can  be  used  to  measure  tilt.  Accelerometer  and  rate  gyro  data  can  be 
combined  using  a  complementary  filter  to  obtain  acceleration  and  tilt  over  the  full  range 
of  frequencies  (Brown  and  Hwang,  1997;  Graber  et.  al.,  2000). 


Waves  can  be  measured  from  the  ASIS  buoy  using  a  six-element  array  of  capacitance 
wave  staffs.  Wave  staffs  are  positioned  on  the  faces  of  each  of  the  five  uprights  as  well 
as  at  the  center  of  the  pentagonal  cage.  An  array  of  this  type  can  be  purchased  from  RBR 
Ltd.  for  $7350.  However,  these  waves  staffs  are  not  robust  enough  for  long  deployments. 
A  new  method  of  making  wave  measurements  from  ASIS  is  currently  under  development 
at  WHOI  and  consists  of  a  similar  array  of  upward-looking  2MHz  Doppler  sonars.  These 
sonars  determine  the  range  to  the  surface  by  emitting  sound  pulses  and  measuring  the 
time  of  flight. 


Item 

Cost  ($) 

Manufacturer 

ASIS  buoy 

35,000 

WHOI 

Surface  mooring  and  tether 

10,000 

WHOI 

D-cell  batteries  (60  kW-h) 

4,000 

Matthews  Assoc. 

SFAS  SODAR 

65,000 

Scintec 

Inertial  motion  unit 

15,000 

BEI  Systran  Donner 

Compass 

1,600 

PNI  Corporation 

Data  acquisition  system 

5,000 

WHOI 

Radar  reflector  &  Marine  light 

1,000 

Carmanah 

ARGOS 

1,000 

Seimac  Corp. 

Freewave  modem  (2) 

2,000 

Freewave  Technologies 

Rain  gauge 

600 

RM  Young 

Temperature  and  RFH  probe 

900 

RM  Young 

Precision  spectral  pyranometer 

2,200 

Eppley  Lab 

Precision  infrared  radiometer 

3,300 

Eppley  Lab 

Doppler  sonar  wave  gauge  (6) 

15,000 

WHOI 

Acoustic  Doppler  Current  Profiler 

25,000 

RD  Instruments 

Tota  1 : 

186,600 

Table  3.3  Typical  ASIS  sensors  and  cost. 
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Chapter  4 

Analysis  of  ASIS  Platform  Dynamics 

4.1  Overview 

In  this  chapter  we  model  ASIS  platform  dynamics  and  analyze  ASIS’s  measured 
response  to  different  sea  states.  We  would  like  to  get  an  idea  of  how  ASIS  responds  in  a 
variety  of  conditions,  so  that  we  can  ultimately  predict  and  understand  how  best  to  deploy 
a  sodar.  To  get  a  good  overview  of  ASIS  response  we  look  at  a  year’s  worth  of  data 
along  the  continental  shelf  from  New  York  and  New  England.  We  use  the  data  to 
perform  a  statistical  analysis  of  buoy  response  to  waves  and  wind.  With  this  data,  we  can 
also  verify  our  ASIS  models  that  were  derived  from  theory.  A  better  understanding  of 
ASIS  platform  dynamics  provides  the  knowledge  required  to  successfully  deploy  a  sodar. 

4.2  Wind  and  wave  conditions 

Wind  and  wave  conditions  used  in  this  analysis  are  chosen  based  on  MVCO  and  National 
Data  Buoy  Center  (NDBC)  data.  A  map  of  the  NDBC  buoys  in  the  Northeast  region  is 
shown  in  Figure  4.1.  The  buoys  selected  to  provide  data  for  our  models  are  shown  in 
Table  4.1. 


Figure  4.1  NDBC  Northeastern  buoys. 
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Buoy 

Latitude 

Longitude 

Location 

44005 

40°1 1'12"  N 

69°09'48"  W 

78nm  East  of  Portsmouth,  ME 

44008 

40°30'00"  N 

69°25'53"  W 

54nm  Southeast  of  Nantucket,  MA 

44017 

40°41'32"  N 

72°02'52"  W 

23nm  Southwest  of  Montauk  Point,  NY 

44018 

41°15'30"  N 

69°17'40"  W 

30nm  East  of  Nantucket,  MA 

44025 

40°15'01"  N 

73°09'59"  W 

33nm  South  of  Islip,  NY 

Table  4.1  Selected  NDBC  buoys  and  their  locations. 


NDBC  buoys  44025  and  44017  are  located  close  to  Long  Island  while  buoys  44008  and 
44018  are  in  deeper  water  southeast  of  Cape  Cod.  Buoy  44005  is  located  near  the  Gulf  of 
Maine.  All  buoys  recorded  wind  speed  and  direction  as  well  as  scalar  wave  properties. 
However,  only  buoy  44025  recorded  wave  direction. 

A  statistical  analysis  of  the  wind  and  waves  at  these  buoys,  and  at  MVCO,  show  that  they 
all  experienced  similar  conditions.  Figures  4.2  and  4.3  present  a  typical  comparison  by 
showing  histograms  of  wind  speed  and  direction,  significant  wave  height  and  mean  wave 
periods  and  direction  (where  available)  for  NDBC  buoys  44025  (off  Long  Island)  and 
44018  (southeast  of  Cape  Cod). 


Therefore,  in  the  following  discussion,  in  order  to  increase  the  sampling  statistics,  we  will 
combine  the  data  from  the  four  NDBC  buoys  listed  above.  Histograms  of  the  combined 
data  set  (from  all  4  buoys)  are  shown  in  Figure  4.4.  Winds  are  seen  to  be  predominately 
from  the  southwest  with  a  mode  of  5  m/s.  Waves  are  predominately  from  the  southeast, 
with  a  mode  of  significant  height  and  peak  period  of  about  1  meter  and  7  seconds, 
respectively. 
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Figure  4.2  Histograms  of  wind  and  wave  conditions  during  2004  by  NDBC  buoy 

44025,  located  off  Long  Island,  NY. 
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located  southeast  of  Cape  Cod,  MA. 


significant  wave  height  (m) 
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wave  direction  (degT) 
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Figure  4.4  Histograms  of  combined  data  from  all  five  NDBC  buoys. 
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From  the  combined  NDBC  data,  we  define  three  sea  state  conditions  as:  “  typical”, 

“high”  and  “extreme”.  “Typical”  conditions  are  taken  to  be  the  mean  wind  speed  of  all 
winds  exceeding  a  threshold  of  3.5  m/s.  “High”  winds  are  defined  as  the  mean  of  the 
upper  5  %  of  winds  above  3.5  m/s.  Wave  periods  and  heights  are  taken  to  be  averages 
over  all  wind  events  in  each  class.  Our  “extreme”  condition  will  be  defined  later  in  terms 
of  a  rare  storm  event.  Table  4.2  shows  a  listing  of  the  characteristics  of  the  “typical”  and 
“high”  classes  defined  from  the  combined  2004  NDBC  data. 


Combined  NDBC  buoy  data 

“Typical” 

“High” 

Wind  speed  (m/s) 

7.6 

14.8 

Dominant  wave  period  (s) 

5.0 

7.9 

Significant  wave  height  (m) 

1.6 

3.3 

Inverse  wave  age,  Ul0  / Cp 

1.0 

1.2 

Table  4.2  Combined  NDBC  buoy  data  for  typical  and  high  classes. 

The  “inverse  wave  age”  is  defined  to  be  the  ratio  of  wind  speed  at  10  meters  to  the  phase 
speed  of  the  waves  at  the  spectral  peak  and  characterizes  the  state  of  development  of  the 
wave  field. 

4.3  Measured  ASIS  response  to  wave  forcing 

In  this  section  we  smooth  and  parameterize  measured  response  functions  in  heave,  pitch, 
roll  and  surge  from  an  ASIS  buoy.  The  response  functions  are  described  in  Graber  et.  al 
(2000)  and  are  for  an  ASIS  similar  to  the  buoy  design  described  in  this  thesis.  In  fact, 
Graber  et.  al.  (2000)  measured  response  functions  for  two  different  ASIS  buoys  and 
showed  that  the  configuration  does  not  significantly  effect  the  result.  The  empirical 
transfer  functions  are  shown  in  Figures  4.5  -  4.8  (E.  Terray,  pers.  comm.).  Each  figure 
shows  the  magnitude  (top  panel)  and  phase  (bottom  panel)  of  the  response  function  (solid 
line),  the  wave  forcing  spectrum  (dashed  line  -  top  panel)  and  the  coherence  (dashed  line 
-  bottom  panel)  between  the  measured  quantity  and  the  wave  forcing.  In  the  case  of 
heave  and  surge,  the  forcing  function  is  taken  to  be  the  wave  height,  whereas  for  buoy  tilt 
it  is  taken  to  be  the  sea  surface  slope.  The  response  functions  shown  are  the  average 
response  over  sea  states  with  significant  wave  heights  ranging  from  0.5  to  3  meters. 
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Significant  Heighl  Range  =  0.5  to  3  m 


Figure  4.5  ASIS  heave  average  response  functions  for  sea  states  in  the  range  of 

Hs  =  0.5-3  meters. 


Significant  Height  Range  =  0.5  to  3  m 


Figure  4.6  ASIS  pitch  average  response  functions  for  sea  states  in  the  range  of 

Hs  =  0.5-3  meters. 
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Significant  Height  Range  =  0.5  to  3  m 


Figure  4.7  ASIS  roll  average  response  functions  for  sea  states  in  the  range  of 

Hs  =  0.5-3  meters. 


Significant  Height  Range  =  0.5  to  3  m 


Figure  4.8  ASIS  surge  average  response  functions  for  sea  states  in  the  range  of  Hs 

0.5  -  3  meters. 


33 


We  parameterized  the  empirical  ASIS  response  functions  by  fitting  each  with  a 
smoothing  spline  (see  Figure  4.9).  For  the  case  of  heave,  pitch  and  surge,  we  continued 
the  transfer  functions  to  low  frequency  using  the  theoretical  ideal  response  to  long  waves. 
So  the  phase  of  the  heave  response  tends  to  zero  and  its  magnitude  is  unity  as  f  0.  The 
phase  of  both  pitch  and  surge  tends  to  jc/2  (the  difference  in  sign  of  the  pitch  is 
definitional)  and  the  magnitude  of  the  surge  tends  to  unity.  We  have  taken  the  magnitude 
of  pitch  to  be  0.5  at  low  frequencies,  in  accord  with  the  observations. 

We  have  included  the  roll  response  for  completeness.  Although  we  expected  the  roll 
response  to  be  closer  to  pitch  than  it  is  (since  ASIS  is  a  regular  pentagon),  it  is  quite 
different.  The  reason  for  this  difference  is  not  known. 
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Figure  4.9  Parameterized  ASIS  response  functions. 
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4.4  Comparison  of  the  modeled  and  measured  ASIS  tilt  response 

In  this  section  we  compare  the  parameterized  ASIS  response  functions  with  data  from  our 
ASIS  deployment.  Since  we  did  not  measure  waves  from  the  buoy,  the  comparison  must 
be  statistical.  We  use  wave  height  spectra  measured  at  the  MVCO  12-m  node  to  estimate 
wave  slope  spectra  for  waves  having  frequencies  below  0.5  Hz.  We  then  combined  these 
with  the  parameterized  response  functions  to  predict  the  buoy  rms  tilt.  Lastly,  we 
compared  the  predicted  buoy  tilt  to  the  rms  tilt  measured  by  the  rate  gyros  on  our  ASIS, 
during  the  same  time  period. 

Ideally,  we  want  to  compare  the  individual  vector  components  of  measured  tilt  with 
predicted  tilt  ( ie .  compare  pitch  to  pitch  and  roll  to  roll).  However,  to  do  this  we  must 
know  the  heading  of  ASIS.  Unfortunately,  due  to  a  compass  malfunction,  heading  data  is 
not  available  from  this  portion  of  the  deployment,  and  so  comparing  the  vector 
components  of  tilt  is  not  possible.  To  get  an  upper  bound  on  the  predicted  pitch  of  ASIS, 
we  multiply  the  wave  slope  spectra  by  the  magnitude  of  the  pitch  response  function. 
Similarly,  to  estimate  an  upper  bound  on  the  predicted  buoy  roll,  we  multiply  the  wave 
slope  spectra  by  the  magnitude  of  the  roll  response  function.  During  this  period  the  wave 
direction  was  approximately  constant  from  the  south.  We  argue  that  because  a  tethered 
ASIS  buoy  makes  a  circle  two  times  a  day  due  to  the  tides,  it  therefore  sees  roughly  equal 
forcing  in  pitch  and  roll.  Hence  the  mean  of  our  two  computations  should  provide  a 
reasonable  estimate  of  tilt.  Figures  4. 10  and  4. 1 1,  respectively  show  the  correlation  of 
predicted  tilt  in  pitch  to  measured  tilt  and  the  correlation  of  predicted  tilt  in  roll  to 
measured  tilt. 
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Comparison  of  Tilts  using  Pitch  Transfer  Function 


Predicted  Tilt  (degrees) 

Figure  4.10  A  lower  bound  on  the  correlation  of  predicted  tilt  to  measured  tilt.  The 
solid  line  is  a  linear  fit  to  the  data  points  and  has  a  slope  of  0.79.  The  dashed  line  is 
1:1.  These  data  have  a  correlation  coefficient  of  0.56. 
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Comparison  of  Tilts  using  Roll  Transfer  Function 


Figure  4.11  An  upper  bound  on  the  correlation  of  predicted  tilt  to  measured  tilt. 
The  solid  line  is  a  linear  fit  to  the  data  points  and  has  a  slope  of  1.3.  The  dashed  line 
denotes  1:1.  These  data  have  a  correlation  coefficient  of  0.52. 

The  mean  slope  of  the  fits  in  Figures  4. 10  and  4. 1 1  is  1 .04,  suggesting  that  the  ASIS  tilt 
response  function  gives  a  reasonably  good  prediction  of  the  buoy  rms  tilt. 

4.5  Modeled  ASIS  spectral  response  to  waves 

In  this  section  we  use  a  model  of  the  wave  height  spectrum  due  to  Donelan  et.  al.  (1985), 
coupled  with  the  ASIS  transfer  functions  given  in  Figure  4.9  to  estimate  the  response  of 
ASIS  in  various  sea  states.  The  latter  are  based  on  the  characterization  of  “typical”  and 
“high”  sea  states  derived  from  NDBC  buoy  observations  on  the  New  England  continental 
shelf  (see  Table  4.2).  We  have  augmented  these  with  an  “extreme”  condition  modeled 
after  a  storm  event  recorded  in  late  2005  at  the  MVCO.  A  summary  of  the  parameters  of 
these  three  conditions  is  shown  in  Table  4.3. 
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Simulation  parameters 

“Typical” 

“High” 

“Extreme” 

U  (m/s) 

7 

15 

25 

Tp  (s) 

5 

8 

10 

Computed  quantities 

U/Cp 

0.9 

1.2 

1.6 

Hs  (m) 

1.2 

3.6 

6.3 

Table  4.3  Wind  and  wave  parameters  used  to  simulate  ASIS  response. 


Note  that  we  have  slightly  changed  the  parameters  defining  “typical”  and  “high” 
conditions  for  better  consistency  with  the  Donelan-Hamilton-Hui  (DHH)  model.  We 
used  our  “typical”  and  “high”  conditions  to  define  parameterized  fetch-limited  DHH 
wave  height  spectra  (Donelan  et.  ah,  1985)  which  we  later  use  for  wave  forcing  in  the 
model.  The  DHH  frequency-direction  spectrum 

F(co,d)  =  ^O(w)/3sec/72{/3[0-  6(co)]}  (4.1) 

where  6  is  the  mean  wave  direction  and  co  is  the  angular  frequency  (rad/s),  and  (3(co) 
parameterizes  the  directional  spread  and  is  given  by 

/3  =  2.6\(co/(op)13 ;  0.56  <  a)/a)p  <  0.95 

/3  =  2.28(co/ft>p)“13;  0.95  <  co/(op  <  1.6 

/3  =  1.24;  otherwise.  (4.2) 

The  non-directional  frequency  spectrum  is 

<f>(u>)  =  ag2w“5(«/wp)expj-(— )4|yr 

a  =  0.006 (Uc  /cpf55-  0.83  <Uc/cp  <5. 

y  =  1.7;  0.83  <Uc/cp  <1; 

7  =  1-7  +  6.0 lg (]UC  /cp);  1  <  Uc  icp  <  5  (4.3) 
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where  r  =  exp|-(co  -  a>p  )2  /2 o2<x>2p  | , 


o  =  0.08[1  +  4/(Uc  /cp)3];  0.83  <Uc/cp  <5  (4.4) 

and  Uc  is  the  average  speed  of  10  m  wind  in  the  direction  of  the  waves  at  the  peak  of  the 
spectrum,  and  cp  is  the  phase  speed  at  the  peak  frequency.  The  parameter  Uclcp,  known 

as  the  “inverse  wave  age”,  is  a  measure  of  wave  development.  The  non-directional  DHH 
spectra  used  to  simulate  our  three  classes  of  wind  and  wave  conditions  are  shown  in 
Figure  4.12. 


Frequency  (Hz) 

Figure  4.12  Wave  forcing  as  defined  by  DHH  spectra  for  “typical”,  “high”  and 

“extreme”  cases. 


The  spectral  response  of  ASIS  to  the  wave  forcing  defined  by  the  DHH  spectra  is  shown 
in  Figure  4.13.  The  rms  responses  of  ASIS  are  listed  in  Table  4.4.  Here,  the  significant 
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response  amplitude  is  defined  as  the  average  of  the  response  to  the  one-third  highest 


waves. 
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Figure  4.13  ASIS  spectral  responses  due  to  DHH  wave  forcing  for  the  three  cases 
(blue  has  values  of  fP  =  0.2  Hz,  U  =  7  m/s;  green  has  values  of  fp  =  0.125  Hz,  U  =  15 
m/s;  and  red  has  values  of  fp  =0.1  Hz,  U  =  25  m/s). 


Significant  Response  Amplitude 

“Typical” 

‘High” 

“Extreme” 

Wave  height  (m) 

1.2 

3.6 

6.3 

ASIS  Heave  (m) 

0.59 

2.5 

4.7 

ASIS  Pitch  (deg) 

0.65 

3.1 

4.6 

ASIS  Roll  (deg) 

0.30 

1.3 

2.3 

ASIS  Surge  (m) 

0.72 

2.7 

5.0 

Table  4.4  ASIS  rms  responses  to  wave  forcing  for  “typical”,  “high”  and  “extreme” 

cases. 
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4.6  Modeled  ASIS  temporal  response  to  waves 

It  is  evident  from  Figure  4.13  and  Table  4.4  that  ASIS  responds  more  to  the  low 
frequency,  energetic,  waves  commonly  seen  in  high  sea  states.  However,  to  understand 
how  well  ASIS  follows  the  surface  it  is  necessary  to  study  the  time-history  of  the  buoy 
response.  To  do  this,  we  use  the  DHH  spectra  defined  in  Figure  4. 12  to  simulate  a 
random  time  series  with  Gaussian  statistics.  This  should  produce  a  good  estimate  of  the 
sea  surface  for  the  “typical”  and  “high”  sea  states.  In  the  “extreme”  case,  non-linearities 
(bound  harmonics)  can  increase  in  skewness  of  the  sea  surface.  Although  in  deep  water 
the  non-linearities  do  not  affect  the  wave  velocity,  the  waves  we  synthesize  in  the 
“extreme”  case  may  slightly  under-estimate  the  forcing  in  heave. 

Figure  4. 14  plots  the  time  series  of  wave  height  in  blue  and  ASIS  heave  response  in  red 
for  the  three  defined  classes  of  sea  state.  As  expected,  ASIS  does  not  follow  the  shorter 
waves  characteristic  of  “typical”  conditions  as  well  as  it  follows  the  longer  period  waves 
that  dominate  the  “high”  and  “extreme”  sea  states. 
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Figure  4.14  ASIS  response  to  time  series  of  sea  surface  height  for  the  three  cases. 

4.7  Wave  overtopping  probabilities 

Our  simulated  time  series  can  also  be  used  to  study  wave  overtopping.  We  simulate 
about  35.6  hours  of  wave  forcing  at  a  2  Hz  sample  rate  for  the  cases  defined  in  Table  4.3. 
In  the  conceptual  ASIS  design  discussed  in  section  3.3,  the  sodar  antenna  box  is  located 
2.3  m  above  the  mean  waterline  (MWL)  of  ASIS.  We  consider  a  wave  exceeding  this 
height  to  be  overtopping.  Our  simulation  produced  no  overtopping  in  the  “typical”  and 
“high”  classes.  However,  in  the  “extreme”  class,  we  predicted  that  eight  out  of  about 
12,800  simulated  waves  would  overtop  the  sodar  antenna  box.  Seven  out  of  the  eight 
waves  were  overtopping  for  a  period  of  0.5  seconds  and  one  was  overtopping  for  a  period 
of  1  second.  A  plot  of  overtopping  waves  for  our  simulated  time  series  is  shown  in 
Figure  4. 15.  The  lower  panel  shows  surface  displacement  relative  to  the  mean  sea 
surface,  with  overtopping  waves  denoted  as  red  stars.  The  upper  panel  shows  surface 
displacement  relative  to  the  MWL  of  ASIS.  Again,  overtopping  waves  are  denoted  with 
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red  stars.  A  dotted  line  marks  the  height  of  the  sodar  above  MWL  (2.3  m)  in  both  panels. 
It  is  noticeable  that  overtopping  waves  were  not  necessarily  the  highest  waves  in  the  time 
series.  This  shows  that  the  phase  of  the  ASIS  response  plays  a  central  role  in  determining 
the  percentage  of  waves  that  will  overtop. 
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Figure  4.15  The  height  of  waves  relative  to  MWL  of  ASIS  (top  panel)  and  the 
corresponding  height  of  waves  relative  to  the  mean  sea  surface  (bottom  panel).  The 
height  of  the  sodar  above  MWL  (2.3  m)  is  drawn  in  with  a  dotted  line  in  both 
panels. 


4.8  ASIS  tilt  response  at  sub-waveband  frequencies 

In  this  section,  we  compute  ASIS  tilts  at  frequencies  below  the  wave  band  to  get  an 
understanding  of  the  buoy’s  low  frequency  response.  At  low  frequencies,  lateral 
accelerations  are  negligible,  so  the  accelerometer  output  at  these  frequencies  measures 
the  actual  tilt. 
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A  model  was  developed  in  order  to  predict  ASIS  buoy  tilt  at  low  frequency  in  response  to 
the  drag  induced  by  winds  and  currents  (see  Appendix  C).  In  order  to  validate  it,  we 
compared  the  predicted  tilt  to  20-minute  averages  of  buoy  tilt  as  measured  by  the 
accelerometer  and  compass  tilt  sensor  during  our  deployment  (Figure  4.16).  Tilts  are 
generally  seen  to  be  small,  less  than  1.5°,  and  agreement  between  the  model  prediction 
and  the  accelerometer-based  estimate  of  tilt  is  within  0.5°.  Differences  between  the 
model  and  measured  tilts  may  be  due  to  the  fact  that  the  wind  and  current  data  used  to 
force  the  buoy  in  the  model  were  taken  from  the  MVCO  12-m  node  and  the  shore  tower, 
whereas  the  ASIS  buoy  was  deployed  ~  3  km  further  offshore  (0.5  km  south  of  the 
MVCO  tower). 


ASIS  Tilting  Model  Comparison 


Figure  4.16  ASIS  sub-waveband  tilt  model  compared  with  data  from  the 
accelerometers  and  the  compass  tilt  sensor. 
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4.9  Wave  band  tilts 

In  this  section  we  use  data  to  examine  ASIS  tilts  in  the  wave  band  for  different  sea  states. 
It  is  important  to  understand  the  distributions  of  measured  ASIS  tilts  in  response  to 
waves.  The  tilts  in  the  wave  band  are  computed  by  high-pass  filtering  and  integrating  the 
rate  gyro  data.  In  this  section,  we  compute  heave,  pitch  and  roll  for  three  sea  states 
characterized  by  the  amount  of  rms  angular  rate:  low  rms  Q  (1.05°/s),  mid  rms  Q 
(1.72%),  and  high  rms  Q  (2.35%).  The  data  used  in  this  section  are  taken  from  our 
ASIS  deployment  (described  in  Appendix  B). 

Figures  4. 17,  4. 18,  and  4. 19  show  histograms  of  measured  ASIS  tilts  in  pitch,  roll  and 
yaw  for  the  three  sea  states  defined  above.  As  expected,  the  distributions  of  tilt  are 
approximately  Gaussian,  with  the  majority  of  tilts  being  under  2°.  A  summary  of  the 
standard  deviations  of  tilt  for  the  three  sea  states  is  given  in  Table  4.5 


Figure  4.17  Histograms  of  pitch,  roll  and  yaw  in  the  waveband  for  the  low  sea  state 


case. 
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tilts  for  mid  sea  state 


0  (degrees) 


Figure  4.18  Histograms  of  pitch,  roll  and  yaw  in  the  waveband  for  the  mid  sea  state 

case. 
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Figure  4.19  Histograms  of  pitch,  roll,  and  yaw  in  the  waveband  for  the  high  sea  state 

case. 


Standard  deviations  of  tilt  (degrees) 

Tilt  direction 

Low  sea  state 

Mid  sea  state 

High  sea  state 

Pitch 

0.56 

0.97 

1.53 

Roll 

0.53 

0.64 

1.05 

Yaw 

0.40 

0.49 

0.70 

Table  4.5  Summary  of  standard  deviations  of  tilt  for  the  three  sea  states  defined  in 

section  4.9. 


4.10  ASIS  drift  response  to  combined  wind  and  wave  forcing 

In  this  section,  we  consider  the  response  of  the  ASIS  “system”  (consisting  of  buoy, 
tether,  and  surface  mooring)  to  combined  wind  and  wave  forcing.  One  tool  with  which  to 
accomplish  this  is  the  ASIS  drift  model  (described  in  Appendix  E).  Unfortunately,  our 
deployed  ASIS  only  carried  an  ARGOS  position  locater  that  could  not  resolve  the  motion 
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of  the  buoy  over  its  50  m  radius  watch  circle.  However,  this  model  can  be  used  to  predict 
the  heading  of  tethered  and  un-tethered  ASIS  buoys  for  given  winds  and  currents.  The 
model  also  predicts  when  the  tether  becomes  slack.  It  is  important  to  know  what 
conditions  will  cause  this  to  happen  because  a  slack  tether  can  lead  to  chaffing  and  a 
premature  separation  of  the  buoy  from  its  mooring  (this  will  be  discussed  in  more  detail 
later  in  this  section). 

Although  we  only  obtained  about  28  hours  of  heading  data  from  our  deployment  of 
ASIS,  due  to  the  malfunction  of  our  PNI  compass,  we  can  compare  these  data  to  our 
model  as  a  check  on  the  model  accuracy.  The  wind  and  current  velocities  used  in  the 
model  were  taken  from  the  MVCO  (12-m  node  and  shore  tower).  The  measured  ASIS 
heading  compared  to  the  ASIS  prediction  for  01-02  October  2005  is  shown  in  Figure 
4.20.  The  data  is  averaged  to  one  minute  and  heading  is  in  degrees  true  in  the  direction 
the  tether  is  streaming  from.  The  vectoral  components  of  the  wind  and  currents  as  well 
as  the  tension  in  the  tether  for  this  time  period  are  shown  in  Figures  4.21  and  4.22, 
respectively.  In  this  case,  the  maximum  tension  in  the  tether  was  approximately  67  lbs. 

It  is  evident  that  the  model  shows  agreement  with  the  data  during  periods  when  the 
tension  in  the  tether  is  relatively  large,  and  when  the  direction  of  winds  and  currents 
align.  The  agreement  is  worse  during  periods  when  winds  and  currents  are  misaligned. 
This  may  be  due  to  rotation  of  the  buoy  about  the  tether  attachment  point  caused  by  the 
misalignment  of  winds  and  currents. 

One  way  to  improve  the  model  in  the  future  would  be  to  equip  ASIS  with  an 
anemometer,  a  downward-looking  ADCP,  and  a  GPS.  These  would  provide 
measurements  of  winds,  currents,  buoy  velocity  over  the  bottom  and  buoy  location,  thus 
improving  knowledge  of  the  wind  and  current  forcing  as  well  as  providing  higher 
resolution  measurements  of  buoy  position  and  motion. 
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ASIS  Heading  Model 


Figure  4.20  ASIS  Drift  Model  comparison  to  data  for  01-02  October  2005  (YD  274- 

275). 
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Tension  (N)  Velocity  (m/s) 


Wind  Vectoral  Components 


Current  Vectoral  Components 


Figure  4.21  Wind  and  current  vectors. 


Magnitude  of  Tension  in  Tether 


Figure  4.22  Magnitude  of  tension  in  the  tether. 
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Since  the  ASIS  drift  model  shows  reasonably  good  agreement  with  the  heading  data,  we 
use  it  further  to  predict  the  ASIS  response  to  combined  wind  and  wave  forcing.  We  are 
interested  answering  two  questions:  [1]  How  fast  does  ASIS  reach  an  equilibrium  speed 
after  a  change  in  forcing?  and  [2]  Under  what  conditions  does  ASIS’s  tether  become 
slack? 

[1]  ASIS  response  to  changes  in  current 

Our  drift  model  is  time-dependent,  so  we  can  use  it  to  compute  the  time  it  takes  ASIS  to 
respond  to  changes  in  the  current.  It  takes  a  free-drifting  ASIS  62.5  seconds  to  reach  a 
velocity  that  is  97  %  of  a  current  velocity  of  0.5  m/s  (current  velocity  is  changed 
instantaneously  at  t  =  0  seconds  from  0  m/s  to  0.5  m/s).  A  plot  of  ASIS’s  velocity  in  time 
is  given  in  Figure  4.23.  A  non-dimensional  time  scale,  x,  for  changes  in  velocity  can  be 
defined  as 

(aU) 

T  =  - - -t 

Mv 

N 

where  t  is  time  (s),  a  is  the  effective  drag  coefficient  ( - j),  U  is  current  velocity 

(rn  /  s) 

(m/s),  and  Mv  is  buoy  virtual  mass  (kg).  A  plot  of  (Vf  -  V)/Vf,  versus  x,  where  V  and 

cl 

Vf  are  the  buoy’s  instantaneous  and  final  velocity,  is  shown  in  Figure  4.24  for - =  1 

Mv 

and  U=  0.1,  0.5  and  1  m/s).  It  is  evident  that  ASIS  accelerates  rapidly  at  first  before 
relaxing  to  its  final  velocity.  The  figure  shows  that  ASIS’s  response  to  currents  collapses 
under  this  scaling.  Since  the  buoy’s  response  time  is  small  compared  to  a  tidal  cycle,  a 
quasi-static  model  is  appropriate  for  predicting  response  to  changes  in  current. 

It  is  important  to  note  that  this  scaling  is  expected  to  break  down  when  the  wind  forcing 
becomes  comparable  to  that  due  to  currents.  Such  dynamical  environments  are  discussed 
in  the  next  section. 
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ASIS  Response  Time  to  0.5  m/s  Current 


Figure  4.23  ASIS  response  to  a  0.5  m/s  current  and  no  winds. 


Figure  4.24  ASIS  response  to  currents  of  0.1,  0.5  and  1  m/s  for  no  winds. 
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[2]  A  SIS  tether  tension 

The  tethered  ASIS  buoy  can  experience  fouling  of  its  mooring  line  when  the  tether 
tension  changes  from  taut  to  slack  and  vice  versa.  This  is  a  potential  problem  because 
fouling  and  chaffing  can  cause  ASIS  to  go  adrift.  Therefore,  it  is  useful  to  compute  tether 
tension  over  typical  a  given  tidal  cycle  for  different  wind  conditions.  We  base  our 
simulation  on  currents  measured  at  the  MVCO  over  a  tidal  cycle  on  1  January  2005.  We 
take  the  wind  direction  to  be  from  the  south,  which  is  typical  for  the  New  England  shelf, 
and  consider  wind  speeds  of  1,  5,  10  and  15  m/s.  The  phase  and  amplitude  of  the 
measured  tidal  velocity  (20  minute  averages)  for  1  January  2005  are  shown  in  Figure 
4.25.  Phases  are  measured  relative  to  the  wind  vector  using  the  usual  mathematical 
convention  (z'.e.  positive  counterclockwise).  The  tidal  velocity  rotates  clockwise  with  a 
peak  velocity  of  0.25-0.3  m/s  and  an  average  magnitude  of  0. 14  m/s. 


T3 


</> 


Figure  4.25  Amplitude  and  phase  of  MVCO  tidal  current  for  1  January  2005.  Phase 
is  referenced  to  the  direction  of  the  wind  vector  (north,  in  this  case). 

The  steady  state  response  of  ASIS  to  the  tide  shown  in  the  figure  above  and  assumed 
southerly  winds  of  1,5,10,  and  15  m/s  is  shown  in  Figure  4.26.  At  a  wind  speed  of  1  m/s 
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the  direction  of  the  tether  mimics  the  tidal  velocity.  As  wind  speed  increases,  the  degree 
to  which  the  tether  follows  decreases.  At  a  wind  speed  of  15  m/s,  the  tether  moves  little 
in  response  to  the  tidal  velocity.  At  this  wind  speed  the  tension  in  the  tether  is  within  2  % 
of  the  drag  due  to  the  wind  alone.  At  winds  speeds  below  5  m/s  the  drag  due  to  current 
contributes  50  -  100  %  of  the  tension.  However,  tension  is  never  that  large.  Even  in  the 
high  wind  speed  case,  the  mean  tension  in  the  tether  is  roughly  100  lbs.  A  summary  of 
the  tether  tension  for  different  wind  speeds  is  given  in  Table  4.6. 


ASIS  Tether  (wind:  1,  5,  10, 15  m/s  toward  North) 


Figure  4.26  Direction  of  ASIS  tether  (top  panel)  and  tether  tension  (bottom  panel) 
for  a  given  tidal  cycle  and  winds  of  1,5,10,15  m/s. 


Wind  speed  (m/s) 

1 

5 

10 

15 

Wind  drag  (lbs) 

0.452 

11.3 

45.2 

102 

Mean  tension  (lbs) 

16.1 

21.0 

49.3 

104 

%  error 

97 

46 

8.3 

1.9 

Table  4.6  Tether  tension  for  different  wind  speeds. 
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Chapter  5 

Sodar  Performance  on  ASIS  and  Conclusions 

5.1  Overview 

As  discussed  in  section  2.4,  sodar  performance  mainly  depends  on  reduction  in  the 
signal-to-noise  ratio  (SNR)  due  to  the  motion-induced  decrease  in  scattering  volume 
overlap  between  the  transmit  and  receive  beams.  Since  typical  SNR  values  over  the 
ocean  are  not  well  known,  we  focus  here  on  the  ratio  SNR/ SNR0,  where  SNR  is  the 
received  SNR,  and  SNR0  is  the  SNR  that  would  have  been  received  in  the  absence  of 
motion.  In  this  chapter,  we  compute  the  probability  of  the  reduction  in  SNR  due  to:  [1] 
heave  (z),  [2]  surge  (x),  and  [3]  pitch  (0),  using  wave  measurements  taken  during  2005  at 
the  MVCO.  We  choose  to  study  these  three  motions  because  they  are  the  generic  building 
blocks  of  more  complex  motions. 

5.2  SNR  probability 

To  analyze  the  performance  of  the  sodar  on  a  buoy,  we  need  the  probability  that  the  SNR 
ratio  exceeds  a  value,  p,  where  0  <  p  <  1,  for  a  given  sea  state.  We  compute  this 
conditional  probability  separately  for  vertical  and  horizontal  translations  (i.e.  “heave”  and 
“surge”)  and  tilt  (i.e.  “pitch”)  as  a  function  of  sea  state.  Averaging  this  conditional 
probability  over  the  probability  of  the  occurrence  of  various  sea  states  gives  the  data 
availability  for  a  specified  reduction  in  SNR. 

To  illustrate  how  the  conditional  probability  is  determined,  we  consider  the  case  of  a  z- 
translation  (“heave”).  The  same  argument  applies  without  change  to  the  other  motions. 
We  are  interested  in  the  probability  density  function  (pdf)  of  a  change  in  buoy  height,  z, 
during  a  time  At.  In  appendix  G,  we  argue  that  the  pdf  is  Gaussian,  and  hence  can  be 
written  as 


P:(z)  = 


-z  I 2(7- 


(5.1) 
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where  the  width,  o,  depends  on  the  sea  state  and  the  time  difference  At.  The  probability 
that  SNR  /  SNR0  >  p  is 


z(p)  / A R  /  \  /  *  d 

PXp  I  a)  -  PZ(SNR/SNR0  >  p)  =  /  ps(z)<fc  /  f  p,(z)dz  =erf(^>)  / erf{ - (5.2) 

o  /  o  CT-V2  /  cr,V2 


where  e/;/(v)  is  the  “error  function”  (note  that  is  has  the  property  erf(y  ->  o°)  ->  1  and 
<?r/'(y  ->  0)  — * *  0),  and  we  have  taken  into  account  the  fact  that  p=0  for  z  outside  the 
interval  |z|  <  A R  (here  A R  is  the  radial  extent  of  the  range  cell).  The  SNR  ratio,  p,  can  be 
determined  as  a  function  of  z  from  an  analysis  of  buoy  motion.  Inverting  this  to  get  z(p) 
permits  evaluation  of  the  probability  P(p|a).  The  form  of  the  probability  in  (5.2)  applies 
as  well  to  surge  and  tilt  ( i.e .  motion  in  x  and  0).  However,  what  must  be  determined  in 
each  case  are  the  functions  z(p),  x(p)  and  0(p),  and  the  widths  cr_,  ox  and  og.  We 
address  this  in  the  next  two  sections. 


5.3  Overlap  percentages 

Although  a  general  motion  of  the  sodar  platform  results  from  the  combination  of  motion 
in  all  six  degrees-of-freedom,  it  is  useful  for  conceptual  purposes  to  consider  how  the 
three  basic  motions  -  heave,  surge  and  pitch  -separately  affect  the  transmit/receive  (T/R) 
volume  overlap  for  a  vertically-oriented  beam.  The  estimates  from  the  analysis  of  this 
simple  geometry  provide  order-of-magnitude  bounds  appropriate  for  both  slanted  beams 
and  more  complex  motions.  We  discuss  these  three  motions  in  turn  below.  In  all  cases 
we  assume  that  the  SNR  ratio  is  equal  to  the  fractional  change,  AR IV ,  in  the  T/R  volume 
overlap,  and  hence  use  these  two  terms  interchangeably. 


[1]  Heave 

If  the  sodar  heaves  an  amount  z  during  the  round-trip  travel  time  (At  =  2 RJc,  where  c  is 
the  sound  speed)  to  a  sample  volume  at  range  R,  then  the  fractional  T/R  volume  overlap 
(shown  in  Figure  5. 1)  is  approximately  given  by 


p(z)  ~  1  - 


A  R 


(5.3) 
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where  A R  is  the  radial  extent  of  the  sample  volume.  This  expression  ignores  the  curvature 
of  the  sample  volume,  and  models  it  as  a  rectangular  disk  of  area  A  and  height  A R.  For 
the  narrow  beams  under  consideration  here  (6dB  fwhm  ~  10°)  this  is  a  reasonable 
approximation.  Equation  5.3  is  inverted  to  get  z(p) 


\z(p)\ 

ct-\/2 


A R 
ctV 2 


(1-P). 


(5.4) 


Substituting  this  expression  into  (5.2)  gives  the  probability  P(p|a). 


1 

0.9 
0.8 
0.7 

o  0.6 

a : 
z 

2  0.5 
a. 

«  0.4 
0.3 
0.2 
0.1 
0 

0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1 

z/AR 

Figure  5.1  Heave  volume  fractional  overlap  (or,  equivalently,  SNR  ratio). 

[2]  Surge 

We  compute  the  fractional  volume  overlap  due  to  surge  by  carrying  out  a  Monte  Carlo 
integration  of  the  overlap  of  the  main  lobe  of  two  translated  beam  patterns.  In  order  to 
suppress  the  effect  of  sidelobes,  we  have  set  the  intensity  distributions  outside  the  main 


Vertical  translation 
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peak  to  zero  ( i.e .  for  values  0  >  5.5°).  Since  the  diameter,  dX,  of  a  range  cell  increases 
with  range,  the  fractional  overlap  in  surge  depends  on  the  radial  distance  to  the  range  cell. 
The  results  of  this  calculation  are  shown  in  Figure  5.2  for  a  number  of  ranges. 

The  curves  in  Figure  5.2  collapse  approximately  to  a  single  curve  when  plotted  versus  the 
“fractional  surge”,  XldX  (Figure  5.3). 


The  dashed  line  in  Figure  5.3  shows  the  model 


X 

~dX 


=  1  -P 


(5.5) 


Although  this  model  slightly  over-predicts  our  Monte  Carlo  results,  the  latter  are  an 
under-prediction  of  the  true  overlap  due  to  the  sharp  cutoff  imposed  to  restrict  the 
calculation  to  the  main  lobe.  As  a  result,  we  will  use  (5.5)  in  our  calculation  of  the 
probability  Px(p  I  crr) . 
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Surge  fractional  volume  overlap 


Figure  5.2  Surge  fractional  volume  overlap  for  different  range  cells  (labeled  next  to 

each  curve). 


Surge  volume  fractional  overlap 


Figure  5.3  Surge  volume  fractional  overlap  versusXMX. 
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[3]  Pitch 

Similarly  to  surge,  the  degradation  in  SNR  due  to  pitch  is  computed  using  a  numerical 
Monte  Carlo  integration.  Two  beam  patterns  are  created,  one  tilted  with  respect  to  the 
other.  A  uniform  distribution  of  random  points  are  selected  in  the  range  0  <  (p  <  2jt, 

0  <  9  <  it  12,  where  0  and  4>  are  the  spherical  polar  and  azimuthal  angles,  respectively. 
The  result  for  tilts  in  a  plane  containing  the  sidelobes  is  shown  in  Figure  5.4. 


Scintec  Beampattern  Overlap 


Figure  5.4  Scintec  SFAS  vertical  beam  pattern  degradation  in  SNR  due  to  pitch  for 

the  shaded  and  unshaded  case. 

Because  this  leads  to  a  multi-valued  function  0(p),  we  make  the  same  simplification  used 
in  the  case  of  surge  by  restricting  0  <  6  <  dc .  As  before,  we  choose  6C  =5.5°  in  order  to 
include  only  the  Monte  Carlo  main  lobe  of  the  unshaded  T/R  intensities.  Figure  2.6 
shows  that  the  main  lobe  is  azimuthally  symmetric,  so  that  this  restricted  calculation 
applies  to  a  tilt  in  any  direction.  However,  due  to  the  sidelobe  structure,  we  expect  that  it 
slightly  under-predicts  the  volume  overlap  for  tilts  in  the  planes  containing  the  sidelobes. 
The  results  are  shown  in  Figure  5.5. 
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Pitch  fractional  volume  overlap 


Figure  5.5  Pitch  volume  fractional  overlap  using  only  the  main  lobe  of  the  unshaded 

beam  pattern. 


5.4  Calculation  of  the  pdf  widths  o 

In  this  section  we  compute  the  standard  deviations,  cr,  of  the  pdfs  of  for  heave,  surge  and 
pitch  increments  for  3-meter  discus  and  ASIS  buoys  as  a  function  of  sea  state  and  time 
lag  At.  We  make  use  of  observations  taken  during  2005  of  the  wave  height  spectrum 
measured  using  the  ADCP  at  the  12  m  MVCO  node. 

A  power  spectrum  was  estimated  for  22,689  20-minute  bursts.  These  were  averaged  into 
20x21  bins  in  mean  frequency  and  significant  wave  height.  Each  averaged  spectrum 
was  used  to  synthesize  a  random  time  series  of  wave  height,  h(t).  Defining  20  uniformly- 
spaced  time  lags,  At,  from  50  ms  to  1  s  (corresponding  to  sodar  ranges  of  8  -  170  m),  we 
computed  time  series  of  heave  increments  dh  =  hit  +  At)  -  hit)  for  each  sea  state.  The 
pdf  of  each  dh  has  the  form  of  a  zero-mean  Gaussian  distribution,  where  a  specifies  the 
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width.  The  overall  the  behavior  of  a  is  found  to  be  very  close  to  linear  in  At  (details  are 
given  in  Appendix  G). 


Spectra  of  wave  slope  were  obtained  from  the  height  spectra  by  multiplying  the  spectral 
value  in  each  frequency  band  by  k2,  where  k(f)  is  the  wavenumber  corresponding  to 
frequency  f  For  computations  involving  the  ASIS,  the  spectra  were  further  multiplied  by 
the  appropriate  ASIS  transfer  function  (described  in  section  4.3).  In  computations 
pertaining  to  the  3-meter  discus  buoy,  we  assumed  that  the  3-meter  discus  buoy  follows 
the  sea  surface  in  heave,  surge  and  pitch  for  all  waves  with  frequencies  below  0.5  Hz,  and 
therefore  used  a  response  function  of  unity  over  the  entire  measured  spectrum. 

Figures  5.6-5.11  show  the  computed  a’s  for  heave,  surge  and  pitch  increments  over 
At=l  second  for  the  3-meter  discus  and  ASIS  buoys.  The  a’s  for  shorter  time  intervals 
would  be  proportionately  smaller.  Sea  states  are  assumed  to  be  characterized  by  the 
significant  wave  height,  Hs,  and  the  mean  frequency,  F p  (computed  as  the  normalized 
first  moment  of  the  wave  height  spectrum). 
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Three  meter  discus  heave 


Fp  (Hz) 


Figure  5.6  The  standard  deviation,  oz ,  of  the  increment  in  heave  over  1  second  for  a 
3-meter  discus  buoy.  Each  curve  shows  the  dependence  of  az  on  the  mean  wave 
frequency  for  a  fixed  interval  of  significant  height.  The  latter  (in  meters)  appears  to 

the  right  of  the  curves. 
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Three  meter  discus  surge 


Figure  5.7  The  standard  deviation,  crx,  of  the  surge  increment  over  At=l  second  for 

a  3-meter  discus  buoy. 
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Figure  5.8  The  standard  deviation,  oe,  of  the  pitch  increment  over  At=l  second  for 

a  3-meter  discus  buoy. 


Three  meter  discus  pitch 
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ASIS  heave 


Figure  5.9  The  standard  deviation,  crz,  of  the  heave  increment  over  At=l  second  for 

the  ASIS  buoy. 
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Figure  5.10  The  standard  deviation,  ox,  of  the  surge  increment  over  At=l  second  for 

the  ASIS  buoy. 
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ASIS  pitch 


Figure  5.11  The  standard  deviation,  og,  of  the  pitch  increment  over  At=l  second  for 

an  ASIS  buoy. 

5.5  Conditional  exceedance  probability 

The  expected  data  availability  from  the  sodar  in  a  particular  sea  state  depends  on  the 
probability  that  the  received  SNR  is  not  too  severely  reduced.  Hence  we  compute  the 
probability  for  SNR/SNR0  to  exceed  a  particular  value  p.  To  do  this  we  use  equation  5.2 
along  with  z(p),x(p),6(p )  (see  section  5.3)  to  obtain  curves  of  the  conditional 
probability  P(p  I  o) .  Plots  of  these  curves  for  vertical  translation  and  tilt  are  shown  in 
Figures  5.12  and  5.13,  respectively.  Figure  5. 14  shows  a  similar  plot  for  the  surge 
calculation,  but  in  this  case  the  curves  represent  o  /  dr ,  where  dr  is  the  diameter  of  the 
main  lobe  of  the  sodar  beam  pattern.  Minimum  and  maximum  values  of  and  o6 

for  both  the  3-meter  discus  and  ASIS  buoys  are  given  in  Table  5.1. 
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min  o  (At=ls) 

max  o  (At=ls) 

min  rms 

max  rms 

3-m  discus  z  (m) 

0.0197 

0.652 

0.0106 

1.10 

x(m) 

0.0195 

0.646 

0.0106 

1.10 

0 (deg) 

0.128 

8.62 

0.151 

5.94 

ASIS  z  (m) 

0.00301 

0.418 

0.00146 

0.887 

x(m) 

0.00621 

0.441 

0.00302 

0.922 

0 (deg) 

0.0629 

3.27 

0.0306 

3.67 

able  5.1  Minimum  and  maximum  a  values  for  the  3-meter  discus  and  ASIS  buoys. 


The  values  of  a  in  Table  5. 1  are  all  computed  for  At=l  second  and  therefore  represent 
upper-bounds  for  all  sodar  range  cells  below  170  m.  Note  that  the  rms  motion  and  a 
values  for  ASIS  are  generally  smaller  than  those  for  the  3-meter  discus. 

The  largest  o,  and  ox  values  in  2005  for  ASIS  were  0.42  and  0.44  m,  respectively, 
corresponding  to  a  significant  wave  height  of  roughly  4  m.  From  Figures  5.12  and  5.14 
we  find  that  for  SNR/SNR0  >0.8  the  effect  of  heave  and  surge  are  likely  to  be  negligible 
for  both  buoy’s  -  even  in  high  sea  states  maximum  values  of  oe  for  both  buoys  are 
considerably  larger  (3.2  and  8.6  for  ASIS  and  the  discus,  respectively).  Then,  for 
example,  if  we  require  SNR/SNR0  >  0.7,  the  data  availability  for  the  discus  is  around 
35%,  where  it  is  70%  for  ASIS.  We  conclude  that  ASIS  performs  much  better  in  high 
sea  states  than  the  3-meter  discus  buoy. 
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Vertical  translation 


0  = 

\4 
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Figure  5.12  The  probability  that  SNR/SNR0  >  p  as  a  function  of  o  for  vertical 
translation,  a’s  are  given  for  a  range  cell  size  i?  =  5m. 


Figure  5.13  The  probability  that  SNR/SNR0  >  p  as  a  function  of  o  for  tilt. 
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Surge 


Figure  5.14  The  probability  that  SNR/SNR0  >  p  as  a  function  of  oldr  for  surge,  dr 

is  the  diameter  of  the  range  cell. 

5.6  Yearly  average  data  availability  due  to  motion 

In  this  section,  we  compute  the  expected  data  availability  over  a  year  for  a  sodar  mounted 
on  both  a  3-meter  discus  and  ASIS  that  is  moving  in  heave,  surge  and  pitch.  We  do  this 
by  averaging  the  conditional  probabilities  P{p  I  a)  for  heave,  pitch  and  surge  over  the 
probability  for  the  observed  occurrence  of  various  sea  states.  Then 

P(P  IA0=  ^Ps(FP,Hs)P(p\<y)  (5.6) 

F„,HS 

where  a  =  o(FP,Hs,\t )  and  Ps  is  the  probability  of  a  sea  state  lying  in  the  bin  centered 
on  ( FP,HS )•  We  compute  a  for  a  lag  of  At  =1  second  since  this  is  an  upper  bound  on  a, 
and  hence  gives  a  lower  bound  on  the  percentage  of  data  availability.  Note,  however, 
that  whereas  P(p  I  o)  depends  on  At  in  the  case  of  pitch  and  heave,  for  surge  it  is 
independent  of  At  (or,  equivalently,  range)  because  the  increase  in  a  with  range  is 
compensated  by  a  corresponding  increase  in  the  diameter  of  the  range  cell.  This  follows 
because  (see  section  5.3) 
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(5.7) 


X(p)  dr 

crtV 2  avV 2 


(1-P) 


where  <ir  is  the  diameter  of  the  range  cell  at  distance  /?.  But  dr~R  and  ox  ~  At  ~  R,  so 
that  their  ratio  is  independent  of  At. 

Plots  of  the  yearly  averaged  data  availability  for  the  three  buoy  motions  seen  by  the  3- 
meter  discus  and  ASIS  buoys  are  shown  in  Figures  5.15-5.17. 


Heave  (5  m  range  gate) 


Figure  5.15  Yearly  average  data  availability  in  heave  for  the  3-meter  discus  and 
ASIS  buoys.  Curves  are  for  the  range  i?=170  m  (Le.  At  =  1  s)  and  hence  give  a 
lower-bound  on  the  percentage  of  data  availability. 
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Figure  5.16  Yearly  average  probability  for  data  availability  in  pitch  for  the  3-meter 
discus  and  ASIS  buoys.  Curves  are  for  the  range  7?=170  m. 


Surge 


P 


Figure  5.17  Yearly  data  availability  in  surge  for  the  3-meter  discus  and  ASIS  buoys. 

Curves  are  independent  of  range. 
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Figures  5.18  and  5.19  summarize  yearly  average  data  availability  in  heave,  pitch  and 
surge  for  the  ASIS  and  3-meter  discus  buoys,  respectively. 


Figure  5.18  ASIS  buoy  %  data  availability  over  2005  for  pitch,  heave  and  surge. 
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3-m  Discus 


Figure  5.19  3-meter  discus  buoy  %  data  availability  over  2005  for  pitch,  heave  and 

surge. 


5.7  Conclusions 

We  have  analyzed  the  operation  of  Doppler  sodar  and  concluded  that  the  principal  effect 
of  platform  motion  is  an  increase  in  the  variance  of  the  Doppler  frequency  estimate  due 
to  a  decrease  in  the  received  signal-to-noise  ratio  (SNR).  The  reduction  in  SNR  is  a 
consequence  of  a  decrease  in  the  spatial  overlap  between  the  transmit  and  receive  beams 
due  to  the  motion.  Based  on  a  full  year  of  observations  of  wave  height  spectra  taken  at 
the  MVCO  12-m  node,  we  have  developed  quantitative  predictions  of  the  probability  for 
the  SNR  (relative  to  the  case  of  no  motion)  to  exceed  any  value  in  the  range  0  to  1 ,  as  a 
function  of  sea  state.  These  conditional  probabilities  were  averaged  over  the  observed 
probability  of  occurrence  of  the  various  sea  states  to  predict  a  yearly  data  availability  for 
a  sodar  deployed  on  either  an  ASIS  or  a  3-meter  discus  buoy  (Figures  5.18  and  5.19). 
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Our  analysis  predicts  that  ASIS  will  be  an  acceptable  platform  for  sodar  measurements  at 
sea.  The  3-meter  discus  buoy  may  also  be  an  acceptable  platform  in  terms  of  yearly  data 
availability,  but  fails  rapidly  to  provide  a  stable  enough  platform  in  high  sea  states. 


SNR  reduction  due  to  heave  and  surge  is  negligible  for  both  the  ASIS  and  3-meter  discus 
buoys.  For  example,  the  yearly  probability  of  data  availability  in  heave  is  80  and  90  % 
for  the  3-meter  and  ASIS  buoys,  respectively.  Data  availability  for  both  buoys  in  heave 
is  similar  to  that  in  surge. 

Sodar  performance  is  most  heavily  affected  by  buoy  tilting  due  to  wave  forcing.  ASIS  is 
substantially  less  reactive  to  forcing  in  pitch  than  the  3-meter  discus  buoy.  For  example, 
ASIS  has  90%  data  availability  for  a  15%  decrease  in  SNR,  while  the  3-meter  discus 
buoy  gives  about  60%  data  availability. 

In  fact,  the  data  availability  models  for  the  3 -meter  discus  buoy  most  likely  over-estimate 
sodar  performance.  In  this  analysis,  the  3-meter  discus  buoy  was  assumed  to  follow  the 
slope  of  the  sea  surface,  when  in  reality,  discus  buoys  are  sensitive  to  mooring  dynamics. 
During  periods  of  strong  forcing  by  winds  and  currents,  their  moorings  can  become  taut 
causing  the  buoy  to  experience  a  more  erratic  motion.  In  contrast,  ASIS  buoys  are 
attached  laterally  to  a  surface  mooring  by  a  buoyant  tether,  and  hence  the  effect  of  the 
mooring  on  the  buoy  motion  is  substantially  reduced.  Also,  because  ASIS  has  a  large 
(~  6m)  draft,  it  can  be  ballasted  to  increase  its  righting  moment,  and  hence  further  reduce 
its  response  in  pitch  and  roll. 

Our  results  strongly  suggest  that  the  combination  of  a  sodar  deployed  on  an  ASIS  buoy  is 
a  feasible  and  practical  technique  for  profiling  winds  over  the  lower  several  hundred 
meters  of  the  marine  atmospheric  boundary  layer. 
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Appendix  A 
Beampatterns 


A.l  Beampattern  theory 

The  Scintec  SFAS  is  a  2  -  dimensional  phased  array  antenna  consisting  of  64  equally 
spaced  sources.  In  order  to  compute  the  beampattern  for  the  SFAS,  we  must  first 
understand  the  case  of  a  simple  1-  dimensional  line  array.  We  want  to  compute  the  far- 
field  intensity  at  an  observation  point  P,  which  is  a  vector  distance  R  from  the  origin.  In 
order  for  this  calculation  to  be  in  the  far-field,  R  must  be  much  greater  than  the  size  of 
the  array.  In  our  array  of  sources,  we  denote  the  nth  source  by  xn .  The  vector  distance 
from  that  source  to  the  observation  point  P  is  denoted  by  rn,  where  rn  =  R  -  x(i.  Each 
individual  source  emits  a  spherical  wave  whose  amplitude  at  P  is  given  by 


exp  (ikrn) 

Pn=an - - 


(A.l) 


where  kn  is  the  wavenumber  of  the  radiation  and  an  is  the  amplitude.  Since  the 
magnitude  of  R  is  much  greater  than  the  size  of  the  array  rn  can  be  approximated  as 


=  |R- 


=  R2- 


2R  »x„  + 


-.1/2 


=  R-  x„  sin  R 


(A.  2) 


so  that 


Pn  =  an 


exp  (ikR) 
R 


exp  (~ikLxn) 


(A.3) 


X 

where  terms  of  order  are  neglected,  and  kL  represents  the  transverse  component  of  the 
wavenumber  k.  The  net  intensity  at  P  is  therefore  given  by 


.  .9. 

A  ' 

-  N 

2 

A 

■ — . 

ii 

ii 

_R_ 

exp  (~ik±xn) 

n=  0 

= 

_R_ 

B(R ) 


(A.  4) 
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The  first  factor  in  (A. 4)  represents  the  fall-off  in  intensity  with  distance  due  to  spherical 
spreading.  The  second  factor 


B(B)  = 


X  rv„exp  (~ikLxn) 

n=  0 


(A.5) 


represents  the  far-field  beampattern  of  intensity,  where  the  weights  {  wn  }  allow  more 
control  over  sidelobe  intensities.  Here,  Ns  represents  the  number  of  transducers  in  the 
array.  It  is  important  to  note  that  {  wn  }  must  sum  to  1  so  that  B(0)  =  1.  Therefore,  the 
far-field  beampattern  is  detennined  by  the  weighted  superposition  of  plane  waves  from 
each  source. 


For  uniform  weighting,  a>n  =  1  /Ns,  the  sum  in  A.5  can  be  computed  analytically.  We 
compute  this  sum  for  an  odd  and  even  number  of  transducers.  For  an  odd  number  of 
transducers,  the  result  is 


B  (#)  = 


s\n(kLNAx /2) 
s\n(k±Ax /2) 


|2 


(A.  6) 


and  for  an  even  number  of  transducers  the  result  is 


B  (#)  = 


2cos((A  +  2)k±Ax  /  4) 


sin(Nk±Ax  /  4) 
sin(k±Ax/2) 


|2 


(A.  7) 


where  kL=  k  sinff .  We  used  the  expressions  in  A.  6  and  A.  7  to  check  our  numerical 
results  (E.  Terray,  personal  communications). 

The  2-dimensional  case  is  computed  similarly  except  the  sources  now  form  a  grid  with 
locations  xn  m  =  (xn,ym)  and  the  observation  point  in  the  far-field,  P,  is  defined  in  terms 
of  polar  and  azimuthal  angles  as 
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R  =  7?(cos(/9sim?,sin(/9sini?,cosi?)  (A.8) 

So  that  (A. 2)  becomes 

rnm  sR-xn  coscpsiniT  - ym  sin  <p  sin#  (A.9) 


and  the  field  at  P  is  now 


Aww. 


expO'Ai?) 

R 


edp(-ik±xn  cos(p)cxp(-ikLym  sin q>) 


(A.  10) 


where  kL  =  Asin# . 


The  intensity  at  P  from  an  Nx  x  Ny  array  is 


/  = 


2 

P 


n= 1  m= 1 


1 

n= 1 


2  wnwm  exp (~ik±xn  coscp)exp(-ikLym  sin<p) 

m= 1 


(A.  11) 


Here,  {kx,k )  =  k±(cos(p,sinq)),  so  B{{k,cp)  can  be  written  in  terms  of  the  1-dimensional 
amplitude  beam  pattern  as 


B(&,<p)  =  \b(kxAx)b(kyAy)\2 


(A.  12) 
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where 


Nx 

b(kxAx)  =  ^wntxp(-ikxxn) 

n= 1 
Ny 

b(kyAy)  =  J  wn  exp (~ikyyn )  (A.  1 3) 

n=  1 


It  is  convenient  to  choose  the  origin  in  the  center  of  the  array,  so  that  for  N  even 


n  =  -N/2  , 


-1,1,...,  N/2 


(A.  14) 


and  for  N  odd 


n  =  -(N  -  l)/2,...,-l,0,l,...(A  - 1)/2 


(A.  15) 


where  xn  =  nAx ,  and  ym  =  mAy . 


A.2  Scintec  Beampatterns 

For  the  Scintec  1-D  and  2-D  vertical  beampatterns  we  use  equations  A.  12  and  A.  13, 
respectively.  We  take  N  =  8  and  Ax, Ay  =  0.5  cm  and  we  multiply  by  the  appropriate 
weights  for  Scintec’ s  shading  pattern. 


The  Scintec  steered  beampatterns  (±19°, ±24°)  are  computed  slightly  differently  than  the 
vertical  beampattern.  The  1-D  steered  beampattern  for  the  unshaded  case  where  f/0  is 
the  angle  of  the  steered  beam  is  computed  using 


b(tf)  = 


5X  exp (~ikxxn  sin(tf)) 

n=  0 


(A.  16) 
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where 


5 


n 


exp {ikxxn  sin(#0) 


(A.  17) 


(Clay  and  Medwin,  1977).  Beampattern  results  for  these  calculations  are  shown  in  Figure 
2.4. 
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Appendix  B 

ASIS/Sodar  Experiment 

B.l  Description  of  experiment 

The  goal  of  this  experiment  was  to  explore  the  use  of  sodar  deployed  on  an  ASIS  buoy  to 
measure  wind  profiles  in  the  MABL.  The  work  was  funded  by  the  Massachusetts 
Technology  Collaborative  in  connection  with  an  offshore  wind  power  initiative,  and  had 
the  long-term  goal  of  enabling  inexpensive  site  assessment  and  monitoring  of  wind 
resources  without  the  need  of  a  meteorological  tower. 

B.2  Instrumentation  and  Data  Acquisition  system 

The  ASIS  buoy  used  in  this  experiment  was  equipped  with  a  variety  of  instrumentation. 

A  Scintec  SFAS  sodar  was  situated  in  the  center  of  ASIS,  on  top  of  a  24-inch  aluminum 
table.  The  sodar  was  equipped  with  a  baffle  that  opened  at  an  angle  of  ~24  degrees  and 
extended  an  additional  two  feet  high.  It  is  important  that  the  sodar  was  the  highest  object 
on  the  buoy  so  that  other  fixed  objects  were  not  visible  in  its  beampattern.  The  sodar  was 
connected  to  the  sodar  processing  unit  (SPU),  which  was  situated  below  the  sodar’ s  table, 
and  was  housed  in  a  watertight  enclosure  box.  In  order  to  measure  buoy  motion,  ASIS 
was  outfitted  with  a  BEI  Systron  Donner  IMU  consisting  of  three  orthogonal  rate  gyros 
and  accelerometers.  A  3-axis  PNI  compass  was  used  to  measure  buoy  heading.  The 
IMU  and  compass  were  located  in  the  same  watertight  enclosure  box  that  housed  the 
buoy’s  Data  Acquisition  System  (DAS). 

The  ASIS  DAS  consisted  of  two  PC- 104  based  stacks  that  recorded  both  analog  and 
serial  data  channels.  One  stack  used  a  Lippert  “Cool  Roadrunner  II”  333  MHz  Pentium 
PC- 104  card  that  was  responsible  for  recording  data  from  all  sensors  except  for  the  sodar. 
Data  was  recorded  on  an  80  GB  hard  drive.  This  stack  also  had  a  Diamond  systems  MM- 
32  analog-to-digital  converter  card,  a  ConnectTech  Extreme  8-port  RS-232  serial  card,  a 
Parvus  4-port  10/100  Mbit  Ethernet  switch,  and  a  Tri-M  IR-104  relay  card.  The  Ethernet 
switch  created  a  local  area  network  within  the  DAS,  and  the  relay  card  provided  a  way  to 
switch  power  on  and  off  to  each  sensor.  A  second  stack  was  necessary  to  log  the  data 
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from  the  sodar.  This  stack  consisted  of  a  Lippert  “Cool  Roadrunner  III”  650  MHz 
Pentium  PC- 104  card,  and  was  connected  to  the  LAN  via  the  Parvus  Ethernet  switch. 
Both  stacks,  along  with  the  compass  and  IMU  were  housed  in  a  Rose  Bopla  watertight 
enclosure  that  sat  on  the  top  deck  of  ASIS.  To  enable  remote  access,  the  DAS  was 
connected  to  a  Freewave  spread  spectrum  900  MHz  modem,  which  was  housed  in  a  Rose 
Bopla  watertight  enclosure.  The  freewave  radio  was  affixed  to  one  of  the  legs  of  the 
sodar  table.  A  picture  of  the  DAS,  IMU,  and  compass  is  shown  in  Figure  B.  1. 

In  order  to  track  the  buoy’s  location,  we  included  both  ARGOS  and  GPS  receivers.  The 
ARGOS  was  an  independent  system  that  ran  on  its  own  batteries.  Both  the  ARGOS 
receiver  and  batteries  are  situated  in  a  watertight  tube  that  was  attached  to  one  of  the  legs 
of  the  sodar  table.  The  Garmin  GPS  receiver  was  attached  to  another  of  the  table’s  legs. 
In  accordance  with  Coast  Guard  regulations,  ASIS  was  also  equipped  with  a  marine  light 
and  a  radar  reflector. 
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Figure  B.l  DAS  and  IMU  enclosure. 


A  layout  of  the  instrumentation  aboard  ASIS  is  shown  in  Figure  B.2.  The  Scintec  SFAS 
sodar  is  situated  in  the  center  of  the  pentagonal  deck  atop  an  aluminum  table.  The 
sodar’s  baffle  is  seated  on  top  of  the  instrument  and  is  supported  by  aluminum  arms  that 
are  welded  to  the  sodar’s  table.  The  SPU  is  visible  below  the  aluminum  table.  The  DAS 
is  the  grey  box  situated  in  front  of  the  SPU.  The  marine  light,  which  is  affixed  to  a  post, 
is  visible  in  the  upper-left  corner  and  the  GPS  (small  white  antenna)  is  visible  in  the 
upper-right  corner.  The  radar  reflector  is  the  larger  white  cylindrical  enclosure  that  is 
partially  obscured  by  the  aluminum  table.  The  ARGOS  is  located  to  the  right  of  the 
DAS. 
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Figure  B.2  ASIS  deck  layout. 


B.2  Power 

ASIS  is  equipped  with  five  battery  cans  located  at  its  foot.  Each  can  houses  490  Alkaline 
D-cell  batteries  that  provide  12  KWh  of  energy  per  can.  This  amounts  to  60  KWh  of 
energy  for  the  ASIS  buoy.  Batteries  are  available  pre-assembled  from  Matthews  Assoc, 
at  ~$800  per  can.  The  power  budget  for  our  experiment  is  about  125  W  (see  Table  B.  1). 
With  60KWh  available,  our  instrumentation  can  be  powered  continuously  for  20  days,  or 
for  longer  at  a  lower  duty  cycle. 

It  is  important  to  note  that  the  sodar  and  its  computer  are  responsible  for  about  67%  of 
the  power  required  for  this  experiment.  This  is  in  part  due  to  the  fact  that  the  sodar 
requires  a  650  MHZ  Pentium  processor  to  collect  data.  Software  modifications  could 
potentially  enable  the  Scintec  SFAS  sodar  to  log  data  in  a  close  to  raw  format,  which 
would  cut  down  on  processor  speed  and  save  about  8  Watts  of  power. 


84 


Power  Budget 

Item 

Manufacturer 

Model 

Power  Req'd  (W) 

333  MHz  Pentium  CPU 

Lip  pert 

Cool  Road  Runner  II 

6.0 

650  MHz  Pentium  CPU 

Lip  pert 

Cool  Road  Runner  III 

15.5 

5  Port  Ethernet  Switch 

Parvus 

Switch  104 

1.5 

20  Output  Relay  Card 

Tri-M  Systems 

IR  104 

0.2 

32  Channel  16  bit  A/D 

Diamond  Systems 

Diamond- MM -32- AT 

2.0 

PC  104  Power  Supply 

Tri-M  Systems 

HE  104 

0.5 

8  Port  Rs-232  Serial  Card 

Connect  Tech 

Xtreme/104  RS-232 

1.0 

Wireless  comms 

Freewave 

Ethernet  radio 

0.3 

DAS  Subtotal 

27.0 

Angular  Rate  Sensors  (3) 

BEI  Systran  Donner 

QRS11 

3.0 

Compass 

Precision  Navigation  Inc. 

TCM-2 

1.0 

Linear  Accelerometers 

Columbia  Research  Labs 

SA-307TX 

0.6 

SFAS 

Scintec 

SFAS 

93 

Instrument  Subtotal 

97.6 

TOTAL 

124.6 

Table  B.l  Fully  instrumented  ASIS  power  budget. 


B.3  Communications 

Communication  with  the  ASIS  buoy  was  enabled  through  the  use  of  two  900  MHz 
Freewave  Spread  Spectrum  Ethernet  Modems.  One  Freewave  was  located  on  ASIS  and 
was  connected  to  the  DAS  via  an  Ethernet  connection.  A  second  Freewave  was  located 
on  the  MVCO  ASIT  tower  and  was  connected  to  a  top-side  node.  An  electro-optic  power 
cable  established  the  communication  of  ASIT  with  the  shore.  The  VNC  (Virtual  Network 
Computing)  software  package  was  used  to  view  and  interact  with  the  computers  on  board 
ASIS  remotely.  This  allowed  for  full  control  of  all  instrumentation  on  the  buoy  from  a 
shore-based  computer.  Should  the  Freewaves  fail,  communication  can  be  established 
with  the  computers  on  ASIS  from  a  ship  via  a  long  Ethernet  cable  that  plugs  into  the 
DAS. 

B.4  ASIS  towing  and  Mooring 

Since  a  fully  instrumented  ASIS  weighs  about  2900  lbs  and  is  10  m  long,  it  typically  is 
deployed  using  a  large-scale  oceanographic  vessel  such  as  the  R/V  Oceanus  or  Knorr. 
However,  the  cost  of  the  use  of  one  of  these  vessels  (>  $19,000/day)  was  not  possible 
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within  our  budget,  so  we  had  to  develop  a  way  to  tow,  deploy  and  recover  ASIS  using  a 
smaller  more  affordable  vessel  like  the  R/V  Tioga,  which  has  a  day-rate  of  $2400.  A 
picture  of  the  Tioga  is  shown  in  Figure  B.3  and  its  specifications  are  listed  in  Table  B.2. 


Figure  B.3  WHOI’s  coastal  research  vessel  R/V  Tioga. 


R/V  Tioga  Specifications 

Dimensions 

Length  60  ft 

Draft  5  ft 

Beam  17  ft 

Space 

Deck  300  ft2 

Lab  200  ft2 

Speed 

Max  -  22.5  kts 

Cruising  -  18.0  kts 

Range  900  mn 

Equipment 

Trawl  winch  (800  m  of 
5/16”  3x19  wire) 

CTD  winch  (750  m 
0.322”  3-cond 
cable) 

A-frame 

Crossbeam:  4600  lbs 
Peak:  2000  lbs 

Table  B.2  R/V  Tioga’s  Specifications. 


86 


Since  the  weight  of  the  buoy  exceeds  the  maximum  weight  that  the  TIOGA’s  A-frame 
can  handle,  it  must  be  towed.  We  have  successfully  towed  ASIS  in  relatively  calm 
conditions  at  9  knots  using  a  line  attached  to  the  buoy’s  foot.  This  towing  line  is  passed 
through  a  sheave  on  the  TIOGA  to  its  trawl  winch.  The  line  is  stopped  off  to  the  deck  of 
the  TIOGA  to  relieve  the  winch  throughout  the  tow.  The  foot  and  batteries  of  ASIS  are 
the  only  things  that  rest  above  water,  so  all  other  instrumentation  must  be  deployed  on 
station.  A  picture  of  ASIS  being  towed  is  shown  in  Figure  B.4.  Once  ASIS  has  been 
deployed,  the  towing  line  is  attached  to  the  main  deck  to  facilitate  recovery. 


Figure  B.4  ASIS  being  towed  at  9  knots. 


Since  ASIS  has  little  reserve  buoyancy,  it  must  be  tethered  to  a  surface  mooring  nearby. 
The  tether  system  used  in  this  experiment  consists  of  50  m  of  Vi  inch  spectra  lined  with 
14  and  6  inch  orange  floats.  The  tether  is  attached  to  a  spherical  surface  mooring.  Since 
the  surface  mooring  is  a  potential  navigational  hazard,  it  too  is  equipped  with  a  marine 
radar  reflector  and  light.  A  picture  of  our  fully  instrumented  ASIS  tethered  to  its  surface 
mooring  is  shown  in  Figure  B.5. 
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Figure  B.5  A  fully  instrumented  ASIS  tethered  to  its  surface  mooring. 
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Appendix  C 
ASIS  Tilting  Model 

We  developed  a  model  to  predict  the  static  tilt  of  the  ASIS  buoy  under  different 
conditions  of  steady  winds  and  currents.  First,  ASIS  was  decomposed  into  several 
components  whose  cross  sectional  area  and  distance  to  the  tether  point  are  determined.  A 
“torque  coefficient”,  D,  for  the  buoy  in  air  and  water  is  then  computed  using  the  formula 

d-\p2c*m.  <C1> 

n 

where  p  is  the  density  of  the  medium  (air  or  water),  Cd  is  the  drag  coefficient,  rn  is  the 
distance  from  the  nth  structural  element  of  ASIS  to  the  tether  point  of  the  buoy,  and  An  is 
the  cross-sectional  area  of  the  n,h  component  as  seen  by  the  air  or  water.  The  torque  on 
the  buoy  in  both  air  and  water  is  computed  from 

T  =  D*  |C/|  *U  (C.2) 

where  T  is  the  vector  torque  of  the  buoy  in  air  or  water,  D  is  the  “torque  coefficient”  in 
air  or  water,  and  U  is  the  velocity  of  the  current  or  wind.  The  torques  due  to  the  air  and 
water-side  flows  are  summed  and  set  equal  to  the  gravitational  torque.  Tilts  are 
computed  via 


6  =  T/(Mg*dr)  (C.3) 

where  0  is  the  tilt  angle  about  the  axis  of  torque,  T  is  the  net  torque  due  to  winds  and 
currents,  M  is  the  mass  of  ASIS,  g  is  the  gravitational  acceleration,  and  dr  is  the  distance 
between  the  center  of  buoyancy  and  the  center  of  mass  of  the  buoy  (CB-CG). 
Verification  of  this  model  is  discussed  in  the  main  text  of  the  thesis  (see  Figure  4. 16). 
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Appendix  D 
Sensor  Calibrations 

The  accelerometers,  rate  gyros  and  compass  were  calibrated  in  the  laboratory  using  a 
mechanically  driven  tilting  table  (Appendix  F).  Each  instrument  was  calibrated  with 
respect  to  its  local  axes,  which  in  turn  are  related  to  a  set  of  axes  relative  to  the  box 
housing  the  data  acquisition  system  (DAS).  The  sensors  were  calibrated  as  assembled  in 
the  box.  A  sketch  of  the  relationship  between  the  various  coordinate  systems  is  shown  in 
Figure  D.  1.  Note  that  the  compass  nominally  has  two  coordinate  systems:  one  for  the 
magnetometers  and  one  for  the  tilt  sensors.  The  magnetometer  coordinates  are  labeled  X, 
Y,  and  Z,  while  the  tilt  sensor  coordinates  are  labeled  B  and  C.  Compass  axes  B  and  C 
are  nominally  oriented  in  the  same  direction  as  axes  X  and  Y,  respectively. 

It  is  important  to  note  that  the  DAS  instrumentation  was  calibrated  as  an  entire  system. 
Data  flow  from  the  sensors  first  passed  through  a  signal  conditioning  board  where  it  was 
low-pass  filtered  (3  dB  roll-off  at  5  Hz),  and  then  digitized  with  a  Diamond  systems  MM- 
32  16  bit  A-D  converter.  So  the  measured  calibration  coefficients  convert  from  the  A-D 
output  (in  bits)  to  the  relevant  physical  units. 
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Forward 


Accelerometers: 

Acceleration  was  measured  using  a  Columbia  Research  Laboratory  SA-307  HPTV  tri- 
axial  accelerometer.  The  accelerometer  was  calibrated  by  a  series  of  static  tilts. 
Accelerometer  axes  B  and  C  were  calibrated  in  the  same  manner,  whereas  axis  A  was 
calibrated  in  a  slightly  different  fashion.  To  calibrate  the  B  and  C  axes  of  the 
accelerometers,  the  box  was  statically  tilted  around  each  axis.  The  tilts  were  measured  by 
a  digital  protractor  (with  a  resolution  of  0.1°)  and  ranged  from  about  -20  to  20  degrees. 
The  accelerometer  output  from  the  digitizer  (A/D)  was  then  averaged  over  the  duration  of 
each  tilt  in  order  to  reduce  noise.  The  average  was  fit  to  the  sine  of  the  angle  of  tilt  using 
a  linear  regression.  The  fits  for  the  B  and  C  axes  are  shown  in  Figures  D.2  and  D.3. 
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Accelerometer(counts) 


2 


Accelerometer  B  Calibration 


4 

x  10 


Figure  D.2  Accelerometer  axis  B  calibration  results  plotted  with  a  linear  fit. 


Accelerometer  C  Calibration 


Figure  D.3  Accelerometer  axis  C  calibration  results  plotted  with  a  linear  fit. 
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We  used  a  slightly  different  procedure  for  calibrating  accelerometer  axis  A.  In  this  case, 
the  DAS  was  statically  rotated  around  axis  C  for  angles  ranging  from  -90  to  90  degrees. 
Tilts  again  were  measured  by  a  digital  protractor.  Since  the  response  for  this  axis  is  a 
cosine,  it  was  necessary  in  this  case  to  rotate  the  box  over  larger  angles  so  that  tilts  in  axis 
A  would  be  large  enough  to  be  seen.  As  before,  the  accelerometer  output  was  averaged 
over  the  duration  of  the  tilt  (which  was  about  four  minutes)  in  order  to  decrease  the 
standard  error  of  the  sample  mean.  The  average  was  related  to  the  cosine  of  the  tilt  angle 
and  was  fit  with  a  linear  equation.  The  result  of  this  calibration  can  be  seen  in  Figure 
D.4. 

.n4  Accelerometer  A  Calibration 


Figure  D.4  Accelerometer  axis  A  calibration  results  plotted  with  a  linear  fit. 

Rate  Gyros: 

Angular  rates  were  measured  using  three  orthogonal  Systron-Donner  QRS-1 1  solid-state 
rate  gyros.  The  three  rate  gyros  were  calibrated  using  our  mechanically  driven  tilting 
table,  which  is  described  in  Appendix  F  of  this  report.  In  order  to  investigate  the 
frequency  dependence  in  the  rate  gyro  output,  axes  A,  B  and  C  were  calibrated  for  three 
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different  periods  of  rotation:  5  seconds,  8  seconds,  and  20  seconds.  Instantaneous 
platform  tilt  was  measured  by  a  potentiometer  (see  Appendix  F  for  details)  attached  to  the 
shaft  of  the  tilting  table.  The  rate  gyros  were  digitized  at  20  Hz  using  the  16-bit  A-D  in 
the  DAS.  We  digitized  the  potentiometer  using  a  separate  10-bit  A-D  that  was 
asynchronous  with  the  DAS  A-D.  Hence  we  could  not  directly  compare  the 
corresponding  time  series,  and  had  instead  to  make  a  statistical  comparison.  A  typical 
rate  gyro  spectrum  used  in  this  calibration  is  shown  in  Figure  D.5. 


Rate  Gyro  Axis  B  Spectrum  for  T  =  5s 


Figure  D.5  Rate  Gyro  Axis  B  spectrum  for  T  =  5s. 

The  conversion  from  A/D  counts  to  rad/s  was  computed  as 


a  = 


/ counts ) 
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where  SQ  {(rad  I  s)2 1  Hz)  is  the  rate  gyro  spectrum,  Se  (rad2 1  Hz)  is  the  potentiometer 
spectrum,  oo  is  frequency  (rad/s),  and  the  integration  is  over  the  dominant  peak  in  the 
spectrum  from  0.15  -  0.2  Hz. 

We  included  only  the  first  harmonic  in  our  calibration  because  digitization  noise  in  the 
potentiometer  spectrum  is  amplified  at  high  frequencies  by  the  factor  of  oo 2 .  This  is 
illustrated  in  Figure  D.6,  which  shows  a  comparison  of  the  rate  gyro  and  potentiometer 
rate  spectra  plotted  on  a  vertical  logarithmic  axis.  We  investigated  the  effect  of  this  noise 
on  the  calibration  coefficients  by  comparing  the  result  to  that  based  on  the  ratio  of  the 
standard  deviations  of  the  differentiated  potentiometer  and  rate  gyro  signals  (e.g.  using 
the  full  spectral  bandwidth).  The  results  from  the  latter  were  inconsistent  with  errors 
ranging  from  10  -  200  %. 


Figure  D.6  A  comparison  of  the  logarithm  of  the  rate  gyro  and  potentiometer  rate 

spectra  vs.  frequency. 
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Conversions  from  counts  (C)  to  angular  rate  (rad/s)  for  the  rate  gyros  are  given  in  Table 
D.l.  Calibration  factors  vary  slightly  over  frequency  by  1-2  %.  The  source  of  this 
variability  is  unknown,  and  we  have  taken  our  calibration  coefficient  as  the  mean  over 
the  three  tilting  frequencies. 


Rate  Gyro  Calibration  Summary 

Axis 

5  Second 
Period 

8  Second 
Period 

20  Second 
Period 

Mean  for  all 
Periods 

Axis  A 

£2  =  2.581x10  "5C 

£2  =  2.505  xlO'5C 

£2  =  2.581x10  "5C 

£2  =  2.556  x  10“5C 

Axis  B 

£2  =  1.547  xlO“5C 

£2  =  1.533  xlO“5C 

£2  =  1.586  x  10“5C 

£2  =  1.554  x  10“5C 

Axis  C 

£2  =  1.598  x  10“5C 

£2  =  1.559  xlO“5C 

£2  =  1.509  x  10“5C 

£2  =  1.555  x  10“5C 

Table  D.2  Rate  Gyro  calibration  summary. 


An  investigation  of  the  rate  gyros  at  rest  was  also  performed  in  the  laboratory.  Data  was 
taken  for  four  minutes  while  the  DAS  was  sitting  flat  on  the  bench.  The  spectrum  of  the 
results  for  each  rate  gyro  axis  is  shown  in  Figure  D.7. 
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Rate  Gyros  at  Rest 


Frequency  (Hz) 


Figure  D.7  Spectra  of  rate  gyros  at  rest. 


A  summary  table  of  the  results  of  the  rate  gyro  data  taken  at  rest  is  shown  below  (Table 
D.2) 


Axis  A 

Axis  B 

Axis  C 

N 

[(rad/s)2/Hz] 

1.353  xl0“7 

8.933  xlO“8 

8.316  xlO“8 

f™ 

[Hz] 

3.928 

5.104 

3.583 

On 

[rad/s] 

1.163  xlO-3 

9.452  x  10“4 

9.119xl0“4 

®  res 

[rad/s] 

8.292  xlO“4 

8.581  xl0“4 

8.526  xl0“4 

Table  D.3  A  summary  of  rate  gyro  results  at  rest. 


where  N  is  the  background  noise  level  of  the  spectrum,  fres  is  the  frequency  of  the 
resonant  peak,  oN  is  the  standard  deviation  of  the  noise,  and  ores  is  the  standard 
deviation  of  the  resonant  peak.  The  peak  is  resolved  by  our  digitization  rate,  but  lies 
above  the  wave  band  of  frequencies,  and  hence  can  be  fdtered. 
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Compass: 

Heading  was  measured  using  a  Precision  Navigation  Inc.  (PNI)  TCM-2  Compass.  This 
unit  measures  three  axes  of  magnetic  field  and  has  a  2-axis  liquid-level  tilt  sensor.  The 
tilt  sensor  axes  of  the  compass  were  calibrated  statically,  by  tilting  the  box  around  each 
axis  for  angles  in  the  range  of  -20  to  20  degrees,  where  the  tilt  angles  were  measured  by  a 
digital  protractor.  The  output  of  the  compass  (which  is  in  degrees)  did  not  fluctuate 
much,  so  data  could  be  recorded  by  looking  at  a  computer  monitor  hooked  up  to  the  DAS 
and  noting  by  hand  the  tilt  angle  recorded  by  the  compass.  The  tilt  angle  measured  by  the 
compass  was  linearly  regressed  against  the  tilt  measured  by  the  digital  protractor. 

Results  for  the  calibration  of  the  compass  tilt  axes  B  and  C  are  given  in  Figures  D.8  and 
D.9. 


Compass  Tilt  Axis  B  Calibration 


Figure  D.8  Compass  Tilt  Axis  B  calibration  results  plotted  with  a  linear  fit. 
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Compass  Tilt  Axis  C  Calibration 


Figure  D.9  Compass  Tilt  Axis  C  calibration  results  plotted  with  a  linear  fit. 

The  SODAR  draws  ~8A  while  running.  It  is  powered  through  two  DC-DC  converters 
located  close  to  the  compass,  and  we  were  concerned  about  the  possible  effect  of  this 
high  current  on  the  measurement  of  heading.  As  a  quick  check,  we  first  looked  at  the 
output  of  the  PNI  compass  while  pointing  it  roughly  North  (as  measured  by  a  handheld 
compass)  inside  the  laboratory.  The  output  of  the  compass  was  recorded  for  several 
minutes  with  the  DC-DC  converters  both  off  and  on  (passing  current  into  a  dummy  8A 
resistive  load).  The  DC-DC  converters  changed  the  heading  by  approximately  10 
degrees  (Table  D.3). 


Compass  > 

orth  Calibration 

State  of  DC-DCs 

Mean  Magnetic  Output 

OFF 

17.452° 

ON 

7.114° 

Table  D.3  Compass  North  Calibration. 
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Because  of  this  result,  we  calibrated  the  magnetic  heading  readout  of  the  compass  over 
360  degrees  with  the  DC-DC  converters  both  off  and  on  (again  with  a  dummy  load).  We 
conducted  the  procedure  outside  in  an  area  having  a  fairly  uniform  magnetic  field.  The 
DAS  was  rotated  in  increments  of  1 5  degrees  on  a  leveled  piece  of  plywood.  Each 
heading  was  recorded  for  several  minutes,  and  the  mean  of  the  instrument’s  heading  was 
compared  to  the  magnetic  heading  measured  by  a  handheld  compass  (resolution  of  1-2°). 
The  calibration  was  performed  with  the  DC-DCs  both  on  and  off.  Results  in  both  cases 
were  slightly  non-linear,  and  therefore  were  fitted  with  a  third-degree  polynomial  of  the 
form 

n=  0 


where  an  are  fitting  coefficients,  H  is  the  handheld  compass  heading  and  /  is  the 
instrument  heading.  The  polynomial  coefficients  for  the  two  cases  are  given  in  Table 
D.4. 


Coefficient 

DC-DCs  ON 

DC-DCs  OFF 

«o 

-2.9358 

-5.2388 

a\ 

0.7315 

0.5249 

a2 

0.0035 

0.0059 

«3 

-1.2887  xl0“5 

-2.0336  x  10“5 

Table  D.4  Calibration  Coefficients  for  the  PNI  compass. 


Plots  of  the  calibration  results  and  the  polynomial  fit  are  shown  in  Figures  D.  10  and 
D.l  1. 
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Compass  Calibration  (DC-DCs  on) 


Figure  D.10  Compass  Calibration  (DC-DCs  on)  data  plotted  with  a  third  degree 

polynomial  fit. 
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Compass  Calibration  (DC-DCs  off) 


Figure  D.ll  Compass  Calibration  (DC-DCs  off)  data  plotted  with  a  third  degree 

polynomial  fit. 

A  summary  of  the  calibrations  of  the  accelerometers,  rate  gyros,  and  compass  tilt  axes  is 
given  in  Table  D.5. 


DAS  Tilt  Sensor  Calibration  Summary 

Instrument  Axis 

Calibration 

Accelerometer  A 

a  =  2.959  xlO“5C  +  0.9964 

tel 

Accelerometer  B 

a  =  2.896  x  10“5C  +  0.0108 

[Si 

Accelerometer  C 

a  =  2.838  xlO“5C  +  0.0097 

tel 

Rate  Gyro  A 

=  2.556  x  10“5C 

[rad/s] 

Rate  Gyro  B 

Q  =  1.558  xl0“5C 

[rad/s] 

Rate  Gyro  C 

=  1.553  x  10“5C 

[rad/s] 

Compass  Tilt  B 

0  =  1.0890,. +  0.1474 

[deg] 

Compass  Tilt  C 

0  =  1.07 10,. -0.6521 

[deg] 

Table  D.5  Summary  of  tilt  sensor  orientations  and  calibrations. 
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Where  a  is  acceleration  in  the  units  of  gravitational  acceleration,  g,  C  is  counts,  9  is 
calibrated  angle  in  degrees,  Q  is  angular  rate  of  change  in  radians  per  second,  and  Ur  is 
the  reported  compass  tilt  angle  in  degrees. 
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Appendix  E 
ASIS  Drift  Model 


The  low  frequency  drift  response  of  ASIS  due  to  the  combined  forcing  of  winds  and 
currents  was  computed  as  follows.  Effective  drag  coefficients  in  air  and  water  for  our 
specific  ASIS  configuration  were  computed  using 

a-\p1c*. A  (E.1) 

n 

where  p  is  the  density  of  the  medium  (air  or  water),  Cd  is  the  drag  coefficient,  and  An  is 

the  cross-sectional  area  of  n'h  structural  element  of  ASIS  in  air  or  water.  The  a’s  were 

N 

found  to  be  0.452  and  657  - j  for  air  and  water,  respectively.  Then  the  air  and  water- 

(: m/s) 

side  drag  on  the  buoy  is  given  by 


D  =  a*|f/,.|*U,.  (E.2) 

where  D  is  the  drag,  and  Ur  is  the  velocity  of  the  air  or  water  relative  to  the  buoy.  We 
also  incorporated  water  drag  on  the  tether  using 

Dr  =  at/,,2  u„  (E.3) 

where  Dr  is  the  drag  on  the  tether,  a  is  computed  as  above  for  the  tether,  Un  is  the 
relative  water  velocity  normal  to  the  tether,  and  u„is  the  unit  vector  normal  to  the  tether. 
Assuming  that  the  tether  extension  is  a  linear  function  of  the  tension,  the  tension  in  the 
tether  is  given  by 


T  =  for  R'>  0 

T  =  0  for  R'<0  (E.4) 
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where  T  is  the  tension,  k  is  the  spring  constant,  R  ’  is  the  tether  extension  beyond  its  un¬ 
stretched  length  ,  and  u,  is  the  unit  vector  tangent  to  the  tether.  The  equation  of  motion 
for  the  buoy  then  is 


M„  —  =  D„  +  Dfl  +  Dr  +  T  (E.5) 

at 

where  D1(,  and  Dflare  the  water  and  air-side  drags,  respectively,  V  is  the  horizontal 
velocity  of  the  buoy  relative  to  the  earth,  and  Mv  is  the  sum  of  the  buoy’s  inertial  and 
added  mass.  It  is  straightforward  to  integrate  equation  (E.5)  for  the  buoy  velocity, 
starting  from  an  initial  value. 
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Appendix  F 
Tilting  Table 

A  mechanically-tilting  table  was  designed  and  constructed  in  order  to  test  the  sodar’s 
sensitivity  to  rotational  motion.  It  also  provided  a  convenient  way  to  calibrate  the  rate 
gyros  inside  the  DAS. 

The  tilting  table  was  designed  by  Ben  Allen  (Woods  Hole  Oceanographic  Institution), 
and  assembled  by  the  author.  The  table  is  made  up  of  two  similar  pieces,  a  base  and  a 
tabletop,  that  are  held  together  by  a  shaft  that  rotates  inside  a  pair  of  bearings.  The 
bearings  allow  the  tabletop  to  move  freely.  The  tabletop’s  movement  is  driven  by  a 
motor  that  is  mounted  on  the  table’s  base.  An  aluminum  arm  connects  the  moving  shaft 
of  the  motor  to  the  top  of  the  table.  The  arm  can  be  adjusted  to  allow  for  a  variable  range 
of  tilt  (the  maximum  is  about  20  degrees).  All  table  parts  are  made  of  aluminum,  with 
the  exception  of  the  top  platform  which  is  made  of  marine  wood.  A  computerized 
drawing  of  the  tilting  table  is  shown  below  in  Figure  F.l. 
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Figure  F.7  A  computerized  drawing  of  the  tilting  table  (figure  courtesy  of  B.  Allen). 

The  motor  is  a  Leeson  DC  straight  shaft  motor  (model  No:  CM34D25NZ61C)  that  is 
driven  by  a  Leeson  DC  control  gearbox.  The  motor  can  rotate  as  fast  as  14  RPM,  which 
allows  for  rotational  periods  in  the  range  of  3  to  20  seconds.  This  is  ideal  for  our 
application  because  this  frequency  range  encompasses  the  band  of  wave-induced  buoy 
motion. 

In  order  to  measure  the  tilt  and  period  of  the  table,  a  Honeywell  RV6  potentiometer  is 
mounted  on  the  table’s  shaft.  A  special  mounting  was  designed  at  WHOI  to  hold  the 
body  of  the  potentiometer  fixed  to  the  base  of  the  table  while  allowing  the  shaft  to  move. 
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The  shaft  of  the  potentiometer  moves  with  the  shaft  of  the  tilting  table.  The 
potentiometer’s  output  is  logged  by  a  10-bit  Dataq  A-D  converter  (DI-194  RS). 

The  potentiometer  was  calibrated  by  statically  tilting  the  table  and  reading  its  output 
using  the  Dataq  A/D  converter.  The  potentiometer’s  voltage  output  divided  by  its 
reference  voltage  was  then  compared  to  the  tilt  of  the  corresponding  calibrated 
accelerometer  as  well  as  the  tilt  measured  by  a  digital  protractor.  Both  calibrations  were 
within  a  few  percent  of  one  another,  so  we  take  the  average  to  be  the  potentiometer 
calibration.  The  result  is 


6  _  _o.595)/ 3.70  x  10 


-3 


V. 


ref 


(F.l) 


where  0  is  angle  of  tilt  (degrees),  Vout  is  output  voltage  in  volts,  and  Vref  is  a  reference 
voltage  in  volts.  Plots  of  the  results  of  the  calibration  using  the  digital  protractor  and  the 
accelerometers  are  shown  in  Figures  F.2  and  F.3,  respectively. 
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Potentiometer  Calibration 


Digital  Protractor  (degrees) 

Figure  F.8  Potentiometer  Calibration  with  the  digital  protractor  plotted  with  a 

linear  fit. 


Potentiometer  Calibration 


Accelerometer  C  (degrees) 

Figure  F.9  Potentiometer  Calibration  with  Accelerometer  C  plotted  with  a  linear  fit. 
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A  plot  of  the  spectrum  of  the  potentiometer  for  a  table  period  of  approximately  8  seconds 
is  shown  in  Figure  F.4.  A  second  harmonic  is  apparent  in  the  neighborhood  of  0.3  Hz. 
This  is  evidence  that  the  table’s  motion  is  not  purely  sinusoidal.  However,  the  ratio  of 
the  area  under  the  peak  of  the  second  harmonic  to  the  area  under  the  peak  of  the  first 
harmonic  is  1.27  x  10“2.  A  similar  second  harmonic  can  also  be  seen  in  the  output  of  the 
rate  gyros  (see  Appendix  D). 

Potentiometer  Spectrum  for  T  =  8  s 


Figure  F.10  A  plot  of  the  potentiometer  spectrum  for  a  table  period  of 

approximately  8  seconds. 
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Appendix  G 

Probability  density  function  of  incremental  motion 

As  discussed  in  section  5.4,  we  are  interested  in  the  effect  on  the  received  signal-to-noise 
ratio  of  three  motions:  heave,  surge  and  pitch.  In  this  appendix,  we  compute  the 
probability  density  function  (pdf)  for  increments  in  heave,  dh,  over  a  time  interval  At 
(; =2R/C ),  corresponding  to  the  round-trip  travel  time  for  a  given  range  cell,  R.  The 
probability  density  for  all  three  motions  has  the  same  form,  and  hence  it  is  sufficient  to 
determine  it  for  heave  alone. 

If  we  know  the  pdf  of  dh  =  h{t  +  At)  -  h(t)  for  a  given  time  lag,  we  can  predict  how  the 
sodar  will  perform  in  different  sea  states.  Unfortunately,  we  did  not  make  wave  height 
measurements  from  ASIS,  and  so  cannot  directly  compute  a  time  series  of  dh.  Instead  we 
use  wave  height  spectra  Sh(f)  computed  from  20  minute  “bursts”  taken  at  the  MVCO  12 
meter  node  continuously  during  2005.  These  spectra  are  used  to  synthesize  time  series  of 
dh  for  each  “burst”,  from  which  we  estimate  a  pdf. 

The  2005  data  set  consists  of  22,689  wave  height  spectra.  In  order  to  keep  track  of 
different  sea  states,  we  group  these  into  20  x  2 1  bins  of  average  frequency,  Fp ,  and 
significant  height,  Hs.  The  grid  is  defined  by  0.088  <Fp<  0.325  (Hz)  and 
0.343  <  Hs  <  4.46  (m)  with  step  sizes  of  0.0124  Hz  and  0.206  m.  Hs  was  taken  to  be  as 
four  times  the  standard  deviation  computed  by  integrating  the  spectral  density  of  the 
wave  height,  while  Fp  was  computed  as  the  first  moment  of  the  height  spectrum.  Figure 
G.  1  shows  a  scatterplot  of  the  22,689  values  of  mean  frequency  and  significant  height 
that  occurred.  We  call  this  a  “Nantucket”  distribution,  because  its  shape  resembles  that 
of  the  island.  Figure  G.2  shows  a  three-dimensional  histogram  of  the  number  of  spectra 
in  each  ( f,Hs )  bin.  We  find  that  168  bins  are  occupied. 
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Figure  G.ll  The  “Nantucket”  distribution  of  28,000  sea  surface  height  spectra. 


22689  20-min  bursts 


Figure  G.12  Pdf  of  sea  surface  height  spectra. 
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For  each  of  the  168  sea  states,  we  synthesize  a  20  minute  random  time  series  (sampled  at 
1  Hz)  of  sea  surface  height,  h(t).  Defining  20  uniformly-spaced  time  lags,  At,  from  50  ms 
to  1  s  (corresponding  to  ranges  of  8  -  170  m),  we  compute  time  series  of  heave 
increments  dh  =  hit  +  At)  -  hit)  for  each  sea  state.  Zero-mean  Gaussian  pdf  of  the  form 

P(dh)  =  — ]_e<dhfna2 

where  o  =  a(Fp,Hs,At)  specifies  the  width,  was  fit  to  each  time  series,  dh.  As  an 
example,  we  show  results  for  the  three  cases  of  high,  mid-range  and  low  sea  states 
defined  in  Table  G.  1  for  the  largest  time  lag  of  At  =  Is.  The  corresponding  average 
spectra  are  shown  in  Figure  G.3. 


Sea  state 

Mean  frequency 
(Hz) 

Significant  height 
(m) 

Time  lag  (s) 

High 

0.10 

4.46 

1.00 

Mid 

0.20 

2.20 

1.00 

Low 

0.27 

1.17 

1.00 

Table  G.l  Summary  of  three  sea  states. 


Figure  G.4  shows  how  close  the  empirical  pdf  is  to  Gaussian,  and  Figure  G.5  shows  the 
standard  deviation  as  a  function  of  time  lag.  These  figures  demonstrate  that  the 
distributions  are  very  close  to  Gaussian  with  varying  widths  (standard  deviations)  that  are 
close  to  linear  in  the  time  lag  At.  Although  we  have  shown  only  3  cases,  the  linear 
dependence  on  At  is  characteristic  of  all  sea  states.  The  three  cases  shown  were  chosen 
because  we  expected  deviations  from  linearity  at  large  values  of  Fp  and/or  Hs.  Small 
departures  from  linearity  were  observed  for  large  values  of  At,  but  overall  the  behavior  of 
a  is  very  close  to  linear  in  At. 
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Figure  G.13  Spectra  of  the  three  sea  states  defined  in  Table  G.l. 


(symbol,  f  ,  H  ,  At):  ('x\ 0.10, 4.4,1),  (’o', 0.20, 2. 2,1),  ('+',0.27,1.2,1) 

P  S 


Figure  G.14  Comparison  of  observations  to  a  normal  (Gaussian)  distribution  for  the 
three  sea  states  and  time  lag  defined  in  Table  G.l. 
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3.5 


Figure  G.15  Standard  deviation  as  a  function  of  time  lag  for  three  defined  sea  states. 


Hence  the  standard  deviation,  a,  of  the  pdf  of  dh  can  be  written  as 


o(Fp,HsAt)  =  /3(Fp,Hs)At . 


(3  was  then  computed  for  each  sea  state  ( Fp,Hs )• 
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